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Communication between peraonnel on the ground aurface and 
people working in underground mine ia very important for normal 
operational work aa well aa for emergency reacue operations , 
Cables are prone to get damaged in the event of a mine disaster. A 
wireless communication system which can transmit and receive 
signals through the earth region over the underground mine area 
plays a crucial role for reacue operations. Normal radio frequency 
electromagnetic (EM) waves get attenuated very rapidly in the 
earth medium, while a very low frequency EM wave can propagate 
sufficient distance with moderate power requirement. This makes 
VLF/ELF ayatema attractive for mine communication. 

For transmitting voice or other signals from the surface of 
the earth to underground mining area, a wire loop antenna may be 
laid on the earth aurface. The loop should be sufficiently large 



to cover a wide nine area. In such a case there may be 
considerable variation o£ current distribution in the loop. 

A receiver operating at VLF/ELF is subjected to atmospheric 
noise which is very much impulsive in nature. Besides, various 
machineries used in mining operation introduce strong interfering 
signals which adversely affect the receiver performance. 

The objective of this thesis la (a) to analyse loop antennas 
on the earth surface, and (b) to study performance of PSK 
receivers in mine environment. A brief description of the work 
follows . 

Chapter 1 gives a brief account of the work in the areas of 
wire antennas for through the earth wireless communication , and 
performance study of receivers In mine environment. 

Chapter 2 describes the derivation of vector (magnetic) and 
scalar potential, due to a horizontal current element radiating 
over the lossy earth. Fourier transform technique has been used 
for obtaining integral representation for the magnetic vector 
potential. The magnetic vector potential has both horizontal and 
vertical components. By applying Larentz gauge to the vector 
potential, the scalar potential for a charge doublet is obtained 
from which the scalar potential for a point charge associated with 
the current element is derived. Both the vector and the scalar 
potentials are represented as Sommerfeld type integrals. The 
horizontal component of vector potential for the special case of 
the source and the observation points being located on the 
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alr/earth Interface haa been obtained in a cloaed form. Thia 
component of the vector potential and the scalar potential due to 
a point charge are the neceaaary potentials for analysing wire 
antennas on the earth surface. 

Chapter 3 describes the formulation of the antenna problem. A 
mixed potential integral equation formulation has been employed. 
The scalar potential does not lend to a closed form representa- 
tion we have used a quasi static approach to obtain a cloaed form 
expression for the scalar potential. The method of moments has 
been used for obtaining current distribution in the loop antenna. 
From the current distribution in the surface loop, the field 
inside the earth haa been computed. 

Chapter 4 describes the derivation of a closed form 
expression for the self impedance of a circular loop antenna on 
the earth surface by considering an uniform current distribution 
for the whole loop. Neumann's formula which is usually used for 
inductance computation haa been suitably modified for the 
computation of the self Impedance. 

In chapter S a method haa been presented for studying perfor- 
iilance of PSK receivers in atmospheric noise which contains an 
impulse-like component superimposed on a homogeneous background. 
This noise haa been modelled as a two component random process. 
One component la Gaussian in nature which represents the 
background noise. The other component is assumed to be a 
stationary Poisson arrival process of pure impulses. 
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A conventional FSK r^eceiver can not perform satlafactorlly in 
presence of strong Interfering signals which may be generated from 
different machineries. The effectiveness of spread spectrum 
receiver for such situation has been studied by considering an 
Interfering tone at the carrier frequency In addition to 
atmospheric noise. 
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CHAPTER 1 


INTRODUCTION 

1.1 THROUGH THE EARTH WIRELESS COMMUNICATION 

Communication to and from miners located deep underground is 
necessary to achieve co-ordinated vork as well as for emergency 
rescue operations. The use of electromagnetic waves for 
underground mine communication, remained almost unexplored till 
late 1960a. The main problem la the earths finite conductivity due 
to which high frequency radio waves cannot propagate through the 
earth beyond a few metres. Using very low frequencies extentlve 
testa were carried out in the seventies at Eagle and Imperial Coal 
Mines in USA. These showed successful voice communication from the 
surface to the mine, and coded Interrupted continuous wave 
transmission from mine to surface [1,2]. 

A downlink (surface to subsurface) transmitter for voice 
communication needs a considerable amount of power, to achieve a 
sufficiently high signal to noise ratio. Since the power 
requirement for surface equipment can be obtained from the mains 
supply, this is not a very serious handicap. On the other hand an 
uplink (subsurface to surface) transmitter should be light, and 
the power requirement has to be met from the miner’s Cap-lamp 
battery. In the event of a disaster, the transmitter should be 
able to continue transmitting signals for sufficiently long time. 
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which can even be two to three days long. To conserve power, 
interrupted continuous wave transmission is usually used. For 
example, every second the transmitter may transmit a pulse for a 
duration of 100 millisecond and remain off for the rest of the 
period. In addition to such "beeps” a set of low bandwidth coded 
messages can also be incorporated in the uplink transmitter. A 
miner receives voice messages from the surface personnel and he 
responds via these coded messages. 

i.&. WIRE ANTENNAS FOR MINE COMMUNICATION 

Both linear and loop wire antennas have been used for voice 
and coded continuous wave transmissions. Linear antenna has been 
experimentally shown to be effective for both downlink [3^4] and 
uplink [5] transmission. The main disadvantage of a linear antenna 
is that it requires some type of grounding at the wire ends to 
provide a contlnous path for the current. Farstad [5] has 
demonstrated the use of roof bolts for grounding of uplink linear 
antennas . 

Ulre loops are easy to implement and can be used both as 
receiving and transmitting antennas. For uplink communication, 
Wait and Spies [6] have suggested the use of a small wire loop 
that can be excited by a portable transmitter. The loop may be 
implemented by laying a wire around a pillar or on the mine floor. 
They obtain the magnetic fields on the surface of the earth due to 
a small hprizontal loop inside the earth by rapreaentlng the loop 
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as a vertical magnetic dipole. Uhen the loop is not very small the 
analysis cannot be carried out by treating it as a vertical 
dipole. Uait and Hill [7] have computed field due to a horizontal 
filamentary current loop in stratified earth, by allowing the loop 
to be of any size, but carrying a constant current. However, it is 
not permissible to assume a constant current in a loop antenna 
when the loop diameter is large. For downlink transmission a loop 
may be laid on the earth surface and it should be sufficiently 
large to cover a wide area inside the mine. In such a loop there 
may be a considerable variation of current distribution. 

In addition to the computation of field, it is also important 
to know the input impedance for the wire loops. Uait [8] has 
computed the self impedance of a rectangular loop on the earth 
surface by considering an uniform current distribution for the 
whole loop. Circular wire loop antennas embedded in a dissipative 
medium have been analysed by King [9]. Comparison of the values 
for input impedance of a small loop obtained assuming constant 
current, and that given by actual analysis, shows that the 
resistive parts differ by an order of magnitude, while the 
reactive parts agree closely. A similar situation is expected when 
the loop is located on a dissipative half space. 

1.3 ELECTROMAGNETIC NOISE IN MINE ENVIRONMENT 

To have a communication link between surface and subsurface 
we have to use ELF or VLF band of frequencies. The performance of 
a communication system at these frequency bands is affected by 
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atmospheric noise which is non-Gaussian in nature. A wide band 
observation o£ ELF/VLF atmospheric noise reveals the presence of 
a non overlapping impulse like component superimposed over a 
homogeneous background noise [10,12]. The impulsive component la 
caused by local thunderstorm activities while the background 
component is due to the numerous unresolved pulses from distant 
world-wide lightning flashes. In addition to the atmospheric 
noise, mach^erles used for normal mining operations create 
Intense electromagnetic interference which is usually dominant at 
the harmonic frequencies of the mains supply. The Interference is 
usually much stronger than the desired signal. 

Two approaches are usually taken to develop an atmospheric 
noise model depending upon the bandwidth of the receiver. At VLF 
the ratio of bandwidth to centre frequency is usually small and 
therefore the received noise may be considered as a narrow band 
process. At ELF this ratio la not small, and so the noise should 
be treated as a wideband process. Omura and Shaft [11] considered 
VLF atmospheric noise as a narrowband process with a log-normal 
envelope. Receiver performance with this noise model was obtained 
by both Intutive and rigorous methods. The rigorous method la very 
difficult to implement. Field and Lewlnstein [12] modelled the 
atmospheric noise as a two component random process- For ELF a 
wide band approach, and for VLF a narrow band approach were taken. 
The model can be used for receivers with envelope detection. 
However, the methodology for evaluation of receiver performance 
with coherent detection as In PSK has not been discussed by the 
authors.,/-/. 
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1.4 PRESENT INVESTIGATION 

Ca3 Loop anLenna on the earth surf^ace 

Uii consider a loop antenna on the earth surface for downlink 
transmission. Fig. 1.1 depicts a downlink communication system. 
The earth is considered to be a homogeneous lossy half space. A 
mixed potential Integral equation formulation has been employed 
for the antenna problem. The method of moments has been used for 
obtaining current distribution in the antenna. The potentials are 
obtained from a horizontal current element radiating over a 
dissipative earth, and they appear as Sommerfeld type Integrals. 
These infinite integrals are highly oscillatory and are difficult 
to evaluate numerically [13,141. When the current element la 
located at the air/earth interface, the Sommerfeld type Integral 
which represent the horizontal component of the magnetic vector 
potential can be evaluated in a closed form. In a mixed potential 
formulation only the horizontal component of magnetic vector 
potential la necessary. The scalar potential is also obtained in a 
closed form, following a quasi-static approach. These potentials 
have been used to obtain current distribution in the loop antenna. 
Magnetic field components inside the earth have been computed for 
the current distribution in the surface loop. Ue also derive a 
closed form expression for self Impedance of a small circular loop 
antenna on the earth surface by considering an uniform current 
distribution for the whole loop. A comparison has been made for 
the input Impedance values obtained by the method of moments and 
those given by this expression. 
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Cb3 A model for atmospheric noise 

We conaider atraoapher ic noiae aa a two component random pro- 
ceaa and follow the wide band approach. One component la the 
background noiae due to diatant and numeroua lightning flaahea. 
Thia component can be deacribed aa a zero-mean Gauaaian proeeaa 
[10,12]. The other component ia due to the local thunderatorm 
activitiea which glvea riae to non-overlapping large apiky 
pulaea. The arrival of theae pulaea ia aaaumed to be a atationary 
Poiaaon proeeaa of pure impulaea. 

Thia noiae model ia uaed to evaluate the performance of a 
coherent PSK receiver. A conventional PSK receiver can not perform 
aatlaf actorl ly in the preaence of atrong interfering aignala , for 
example the harmonic int erf erencea in the mine environment. To 
effectively auppreaa the interfering affect, apread apectrum 
communication technique can be employed [15,16,17]. The perfor- 
mance of a direct aequence apread apectrum receiver in the 
preaence of a tone interference at the carrier frequency beaidea 
the atmoapheric noiae, haa alao been atudied in this Inveatlgatlon . 



CHAPTER 2 


HORIZONTAL CURRENT ELEMENT OVER DISSIPATIVE EARTH 

a.l INTRODUCTION 

This chapter deacrlbea the derivation o£ vector (raagnetic) 
and acalar potential for a time harmonic electric dipole source 
located over a dissipative half space. Fourier transform 
technique [18,19,20] is employed for obtaining an integral 
representation for the magnetic vector potential. From the 
magnetic vector potential, the scalar potential of a point charge 
associated with the horizontal dipole is obtained. These 
potentials are used in Chapters 3 and 4 for the analysis of a loop 
antenna laid on the surface of the earth. 

The effect of finitely conducting earth on the electro- 
magnetic field excited by a time harmonic dipole was first 
analysed by Sommerfeld [21]. Banos [22] follows an unified 
approach irrespective of the dipole being electric or magnetic, 
vertical or horizontal, located in the air region or embeded in 
the conducting medium. Special features such as nonlsotroplc and 
stratified medium have been discussed in the literature [23, 24]. 
Hertz vector potential has traditionally been used to express the 
fields [22,23,24]. Ue express the electric field in terms of the 
magnetic vector potential and the scalar potential due to a point 
charge associated with the current element. By using the vector 
and the scalar potentials, the electric field due to the current 
element can be expressed in a two component from, each component 
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being identified separately with the current and the charge. In 
Chapter 3 we analyse a loop antenna on the earth surface by using 
a mixed-potential (vector and scalar) integral equation 
formulation. In this Integral equation only horizontal component 
of the vector potential la necessary, although the vector 
potential has both horizontal and the vertical components. It will 
be shown In the present chapter, that the horizontal component can 
be expressed in a closed form. In Chapter 3 the scalar potential 
la expressed in a closed form by using a quaai-atatlc approach. 
The quaal-atatic approximation la defered to Chapter 3, since it 
depends on loop dimension. In this chapter we consider only an 
isolated current element. 

2.2 HELMHOLTZ EOUATIONS 

A time harmonic horizontal current element over a dissipative 
half apace la shown in Fig. 2.1. The dipole la oriented in the 
direction of positive x axis and is located at a height z^. The 
source is restricted to be in the air region only, and in the 
limiting eaae, z^ is allowed to be zero when the dipole is located 
at the interface. For a horizontal electric dipole, the boundary 
conditions cannot be satiafied with only the x component of 
magnetic vector potential unless the dipole la in unbounded apace 
or above a perfectly conducting surface. It is necessary to 
include a z component of the vector potential to satisfy the 
boundary conditions. For a x directed current element of dipole 
moment pCp=IdT), with a delta function distribution one can write 
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Fig. 2.1 A current element over lossy earth. 
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Time dependence of e^“^ ia asaumed. The position vectors R and r; 
represent the field point co-ordinate Cx,y, 2 ) and source point 
co-ordinate respectively. The propagation constants 

and are given by 
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where permeabilities and of the two regions have been 

assumed to be same, and {i = 0,1}, and a> denote the 

permltlvltles and angular frequency respectively. The permitivity 
of the lossy half apace can be written as 
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where is the relative complex effective dielectric constant and 

can be written in terms of relative dielectric constant and 

r 

conductivity <y as 



The field components in terms of the vector potentials are 
given by the following relations. 
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where for air £ = £^, and for the lower half space 
The following Fourier transform pair is defined. 
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where £ and f represent any component of the vector potential in 
the space and the transform domain respectively. Application of 
Fourier transform to the set of Helmholtz equations given by (2.1) 
leads to a new set of Helmholtz equations, but having only one 
Independent variable, namely the z co-ordinate. They ^re as 
follows . 
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with Re(>'^) S: 0 and Iiii(ik^) ::S 0 , (i = 0,1). 

The solutions for mai^netlc vector potentials in transform 
domain can be written as 
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a. 3 BOUNDARY CONDITIONS 

The boundary conditione at the interface (2=0) in the apace 
aa well aa in the tranaforoi dotnaina are aa followa. 

At z - O : 


(a) Ey. ia continuoua, hence from (2.5e) one e^ta 
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which can be written in tranaform domain aa 
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(b) ia continuous, hence from (2.9a) and (2.5d) one can write 
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(c) Continuity o£ and (2.5a) glvea 
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(d) From continuity o£ H and (2.5b) one gets 
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Besides these boundary conditions at the Interface 

following conditions are to be satisfied at the 

location z = z . 
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(e) is continuous, and hence from (2.9e) one can write 
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(£) Integrating left and right hand aidea of (2.7a) from to 
one geta 
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2.4 IMPOSITION OF BOUNDARY CONDITIONS 

By impoaing the boundary conditiona, a aet of equationa for 
the unknowna A, B, C, D, E and F can be obtained aa followa. 


At z B O : 
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B + C = E . 


(2.10a) 
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At z ■ z : 
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From (2. 9k), (2. 8a) and (2. 8b) 
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From the equations (2.10a), (2.10b), (2.10d), (2.10e), (2.10£) and 
(2.10e), A, B, C, D, E and F can be obtained. They are given as 
follows . 
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F == D 


(2.11f ) 


where R = 


^Q~^l 


C2.ilfi) 


and 


M = e 




C2.11h) 


2. S MAGNETIC VECTOR POTENTIALS IN SPACE DOMAIN 

Taking Fourier inverse transform of the set of equations 
C2.Sa) to (2.8d), the vector potentials in space domain are 
obtained. They can be written as 
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In (2.12a) A has been written In the form 
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where R^ = (x-x^)^ + (y-y^)^ + (z-z^)^ 
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one can write A and A in the following form 
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where rJ = (x-x^)^ + + (z+z^)^. 
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1 IT * 

and can be identified aa the direct and the image fielda 

when the dipole ia located over a perfectly conducting plane. 
can be identified aa the neeeaaary correction term for an 
imperfect ground. Combining A^ , and A*^ and aubatltutlng D, E 

and F from equatlona (2. lid), (2.116) and (2. Ilf), the magnetic 
vector potentlala l.e., equations (2. 12a) to (2.12d) can be 
rewritten aa follows : 


A Cx,y,z) 
xo 




-Jk R 
•' o 




-Jk R. 
■' o 1 


4nR 


4nR. 


rr" W 

Jj 

-00 


(z +2 ) 


z ^ 0 


doij^d«2 > 


(2. 14a) 


00 


A^l(x,y,z) 


ff- 
j4n^ JJ 
■' -00 


N 6 d«.dtt. 

1 2 


+y. 


z 0 (2.14b) 


^ jj 

“*00 


- (z+z ) 
N»t o o' 

PP„ rr ) N 6 


i 2 ^ i Z 

Vo ■" Vi 


“l 


z i 0 


(2.14e) 


J(y.z-y z ) 
00 . •'^'^1 *^ 00 ' 

rr C>'i-:r«)N e 


PPo rr ^^1 " 

(x,y,0 = — 2 JJ -2 —2 

j4n JJ k. y + ky. 

-00 1' o o 1 


where N is given by (2.13e) 


1 


z < 0 


( 2 . 14<|) 


The integrals in the expressions (2.14) to (2.17) are in the 
form of surface integrals over the entire ^2^ plane. 
These integrals are basically of the same form, excpet (2.16) and 
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(2.17) which contain the additional term The term can be 

removed by defining the following relatione. 


A = U 
zo dx zo 

A . = u . 

zl dx zl 


(2,15) 


where W and U . do not contain at. and are given by 
zo zl 1 


U.. =— 4JJ , - 3 . C2.16) 

-00 1 


zo 


An 


y + k y. 
1^0 o*^ 1 


W 


zl 


-PP 


4fi 


00 






jC^.z-y z ) 


-00 


1,2 ^ .2 
* Vl 




(2.17) 


The surface integrals can be expressed in form of cylindrical 
co-ordinates both in the space and the transform domains using the 
relations 


(x-x^) = r coa0 
CWq) = ^ ain0 
r = (x-x^) + (y-y^) 


for the space co-ordinates x and y, and 
= X cos/9 


^ ain/9 


.2 2^2 
X = «, + 

1 2 


for the transform variables and 


(2. 18a) 


(2.18b) 
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From (2.76) and (2. 18b), y and y. can be written as 

o 1 



with ReC?'^) & 0 and ImCy*^) <0, {i = 0,1). 

The area doij^doi^ can be written as XdXd/9 in polar co-ordinates. 
With these transformations N can be written as 


j(a + o^Cy-y,,)) 

N = e 

j{Xrcoa0 coaft + XrsinC^ sin/?) 

as e 

jXrcoa(/?-0) 



= e 



• 


(2.20) 

From (2.16) 

and (2.20) 

one 

can write 



U (r , z ) = 
zo' ’ 


' (r 

l'^o> 

- iy (z+z ) 
o' o' 

e 

2n 

XdX 1 e 

0 

jXrcoa(/?-0) ^ 

4n2 J 

0 









(2.21) 

Using the 

identity 

for 

the Bessel 

function 

of zero order 

[19], 







Jo(Xr) 

1 

ft 

Zn 

t • 

jXrcoa(/?-0) ^ 


(2.22) 
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one can write 


W (c,z) 
zo^ * ' 


-jy Cz + z ) 

to, . o 

"P/J- r e 


2n 


M ^ 


1,2 ^ , 2 
k.y + k V. 
o 1 


J^(X.r)XdX 


(2.23) 


can be obtained from (2.23) by taking the x derivative as 
defined by (2.15). Doing similar exercises one obtains the 
magnetic vector potentials as follows. 


A (r , z) 
xo^ ’ 


PP^e 


Jk^R 


PPo® 


Jk^R, 


00 


4nR 


- jy (z+z ) 
pfj p o o ^ 

x-D + 1 J (Xr)XdX , 

4nR^ j2n J 


z > 0 


(2.24) 


A (r,z) = 
zo^ * ' 


pp 

2n 


dx J 


-jy (z+z ) 
00. ^ ■'^o o"^ 


0 


J (Xr) XdX, z > 0 
o 


(2.25) 




pp fT ■'^^1“ -^o-o^ 

A ,(r,z) = T J (Ar) XdX, 

xl^ ' J2n J ^ 


z < 0 


(2.26) 


A^lCr.z) = 


'=f'o 5_r 

2n dxj 


00 




i(r^z-y^z^) 


1,2 ^ ^2 
k.y + k y, 
o o' 1 


J^(Xr) XdX. , z < 0. 


(2.27) 


These integrals, which are of the Sommer f eld type, in general, can 
not be evaluated in a closed form. 
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a. 6 CLOSHED FORM EXPRESSION FOR A Cz » z a 05 

xo o 


Uheti the source and the observation points are located at the 
interface the horizontal component of the vector potential has a 
closed form expression which can be obtained as follows. From 
(2.24) and (2.26) when z = z^ = 0, one can write 

A (r,o) = A ^(r.o) 
xo^ ’ xl ^ 


j 2n 


QO 


J (Xr)XdX 


I -JEL 

0 ° ^ 


(2.28) 


Multiplying the numerator and denominator of (2.28) by 


can write (noting that y , = r * ^ 


- X**) 


A (r.o) 
xo ' ' 


^^o r 
T2rF J 


00 


Cy^-y^) J^(kr)XdX 


0 


1 2 .2 

k. - k 
1 o 


(2.29) 


-T-T- {"^4 

2r(k‘-k‘) L J 


00 

where 11= j J (kr) XdX (2.30) 

J o o 
0 

00 

V = -■ J ■ (2.31) 

■*0 

The integrals (2.30) and (2.31) are of the similar form and hence 
only one of them will be treated for closed form derivation. 
Differentiating (2.30) w.r.t, k^ one obtains 
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OD 


<Xl 


k J iXv^XdX 


j r p K. U (.A. 

k- J XdX = J f 

o J „ J ^ ' o 


(2.32) 


The r.h.a of (2.32) has a closed form expression [19] which gives 


CfCl 


^r— Jr J (Xr) XdX = 
dk o o^ 

o J „ 


k e 
o 


-jk r 
o 


(2.33) 


Integration of (2.33) w.r.t. k^ gives 


OCr 


-Jk r 


J 


Jy J (Xc) XdX ~ 


p k e 


— dk 
r o 


-jk r . 

e o .Is 

2 — ^-^^o r^ 


(2.34) 


Since, 


-Jk r 

i ) 

r dr r 


- Jk r 


(2.35) 


one can write 


00 


J 


(Xr) XdX = ^ 


Id .e 

5? 


-jk r 
o 


(2.36) 


Doing similar exercises one gets from (2.31) 


00 


J -J^'iJ^(Xr) XdX = 


- jk.r 

i ^L. (.f ) 

r dr ' r 


(2.37) 
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Hence from ( 2 . 36 ), ( 2 . 37 ), ( 2 . 30 ), (2.31) and (2.190) one can write 

the final closed form expression for A as 

xo 


PM. 


A (r,o) 
xo ’ 




1 £_ 

r dr 


-Jk r 
6 o 


~ Jk. r 
e 1 


(2.38) 


2.7 SCALAR POTENTIAL DUE TO A POINT CHARGE 

From the magnetic vector potential A^, a scalar potential 

can be obtained by using Lorentz gauge. The potential is the 

scalar potential associated with the dipole (current element). The 

+ 

scalar potential due to a point charge q associated with the 
current element can be obtained from by the relation [26] 

^ ^ — (^*) ( 2 . 39 ) 

^O dxo 


where the derivative is to be taken at the source point along the 
direction of the current element. Using Lorentz gauge one can 
write 



(2.40) 


For the horizontal current element the vector potential A^ has x 
and z components , and hence from (2.39) and (2.40) one can write 





A + 


xo 



(2.41) 
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By integrating (2.41) w.r.t. source variable one obtains 



11 

+ o 

2 j + 

4> + 

X z 


(2.42) 

where 


^ J fl; 

xo 

(2.43) 


*2 = 

71 f <li 

xo 

(2.44) 



xC 

o 



Using (2 

+ 

.43) and (2.13a), can be written as follows ; 




+ + + 
^xl ‘^x2 


(2.45) 

where 

Ki 

^2 J ^^x xo-* 
o 

dXo , 

(2.46) 


^x2 

2 J ^ax 

dx^ . 

(2.47) 



MS. 

o 




4,* 

^x3 

= > r f£_ a'^ ) 

^2 J '^^x xo'' 

dx 

o 

(2.48) 


Using (2.13b) one can write r.h.s. of (2.46) as 




^ A* ) dx = 
dx xo*^ o 


> r 
J 


a 

6>x 




oa - 
e 


}6ti‘ 


-on 


>qI(z-Zo)| j^„^(x-x^)+oi2(y-y^)} 


dotj^doi2 


dx 


(2.49) 
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The derivative w.r.t. x introduces a term ja^ while the integra- 
tion w.r.t. X introduces the term — r~- This gives 
o 


> \c^ A 

k2 J'^dx ^ 


) dx 
xo”^ o 


o 


xo 


Similarly one gets 


o 




i A ) dx 
dx xo-^ o 


a"" 

2 *o 
o 


r 

J 


c|- ) dx 

^dx xo o 


■is A = 

1^2 xo 
o 


(2.50) 


(2.51) 


(2.52) 


Using the set of equations (2.45) to (2.48), and (2.50) to (2.52) 
one can write 



(2.53) 


Substituting for A from (2.24), one can write (2.53) In the form 

xo 


X 


qdle 

4nR& 


Jk R 
o 


qdle 


-JVl 


4nR7c 
1 o 


qdl 

line 




oe -j^'^(z+z^) 

— J (Xr)XdX, 

0 ° ^ 


(2.54) 


where use of the continuity equation I = ji.>q and the relation 
1 

— = — = — have been made, 
k. o 
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From C2.l4c)# the derivative -s— A can be written as 

&z zo 




rr 

J 4 n^ JJ 

-00 




1,2 2 


a^da^da^ 


(2.55) 


where N is given by (2.136). Integration of (2.55) w.r.t. gives 


-jy (z+z ) 

® , ... ■'o' o 

PM,, rr 


J^^z " . 2JJ .2 ..2 

J j 4 n 


da^da2- (2.56) 


The double Integral in (2.56) can be transformed Into cylindrical 
co-ordinates as described in section (2.5). Multiplying (2.56) by 


j6i 

o 


and comparing with (2.44) one gets 


-y (z+z ) 

.. CO . . 'o' o' 

ddl r J'l ) ® 

- - T - rs 

k,|K + k y* 

1' o o' 1 


- - W 

Qcll p 
]2n£ 


(2.57) 


Substituting for and ip* from equations (2.54) and (2,57), 

X 2S 

respectively, one gets 


4^ 


-jk R -jk R 

qdl 6 ° qdl e 


OD 


4nR<? 


where 


4nRi^^ 




qdL.r ^ 


® ^ 0*^1 k“y^ + k~yj^ 


~ \y (z+z ) 

e ° ° J^(Ar)AdX 


( 2 . 58 ) 

( 2 . 59 ) 
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Uaing the relation y 


i ' tAT ^ 


X , w can be written as 
o 




- kj 

o 1 


1 . - *'1^ 

,2 ^ , 2 
k.y + k 
o o 1 




k ) 

0^0 * 1 ^ 


.,2 , 2 . 2 ^ 4 2 . 
(•'o - ‘'1> cVo * Vi> 


k V + k y . 
1 o o' 1 


(2.60) 


Substituting for from (2.60) in (2. 58) we get 




- jk R -jk R. „ s /■ . s 

.. ■' o .. o 1 .,2 00 - iy Cz+z ) 

qdl e qdl e qdl<o /J r o'^ 


4rrRs 


AnR^ j£? 

1 o 


j2n 


I r o 

O I 6 

“ J. 


k y +k y 
0 1 o o 1 


^ J (Xr)XdX. 

M O 


(2.61) 


which is an exact expression for the scalar potential In 
Chapter 3, this expression is used to arrive at a closed form 
expression for the Green’s function for a point charge associated 
with a loop antenna located on the earth surface by employing a 
quaslstatic approximation. 



CHAPTER 3 


LOOP ANTENNA ON THE EARTH SURFACE 
FOR UNDERGROUND MINE COMMUNICATION 

3. 1 INTRODUCTION 

High frequency ulrelesft communication through the earth is 
not feasible due to strong attenuation of electromagnetic fields 
while propagating through the conducting earth. By using base-band 
modulation voice communication can be achieved for downlink, l.e., 
surface to subsurface communication [2]. A downlink transmitter 
needs a considerable amount of power, which can be provided by the 
mains supply. A loop antenna may be laid on the surface of the 
earth to provide a communication link from surface to underground 
mine area. The loop should be sufficiently large to cover a wide 
area inside the mine. In such a loop there may be considerable 
variation of current distribution along the loop periphery. 

To analyse wire antenna over the dissipative earth both 
rigorous [27 J and approximate techniques [28,29,30,31] have been 
employed. The approximate techniques break down when the distance 
between the antenna and the earth surface is much smaller than the 
wavelength at the operating frequency. The rigorous method is 
based on Sommerf eld-integral formulation, however, the integral is 
to be evaluated numerically and it is very time consuming. In the 
present investigation the problem of a loop antenna on the earth 
surface (Fig. 3.1) has been formulated in terms of Sommer f eld-type 
integrals. The analysis is carried out by the method of momenta. A 
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A mixed potential integral equation formulation has been employed 
[27,32]. For the magnetic vector potential we shall use the closed 
form expression given by (2.38). The scalar potential la also 
obtained in a closed form following a quasi-static approach. These 
potentials are used to obtain current distribution in the loop 
antenna. Magnetic field components inside the earth have been 
computed due to the current distribution in the surface loop. 

3.2 FORMULATION OF THE ANTENNA PROBLEM 

To analyse a horizontal wire antenna over an Imperfectly 
conducting ground, first the formulation of the problem is carried 
out assuming the antenna to be in free space. Necessary 
modifications are then incorporated in this formulation to allow 
the presence of the dissipative half apace. Fig. 3.2 depicts a 
horizontal antenna of arbitrary shape situated at a height h from 
the ground surface. The wire la assumed to be a perfect conductor. 

3.2.1. Wire antenna In free space 

For an incident electric field E. , the scattered field 

1 a 

from the conductor can be written as 

1 = -JtoA - (3.1) 

3 

where, the magnetic vector potential A and the scalar potential ip 
are due to the current and charge distributions on the surface of 
the conductor. Since, the wire radius is much smaller than the 

wavelength, thin wire approximations can be made according to 




Fig.3.1 Circular wire loop on the surface of the earl 
(a) Top view 
( b] Sectional view at X-X 


Y 



Fig.3.2 A horizontal wire 

(a ) Plan view on 

(b) Side view at 


t 

K-U -1 




air b 


y jT ^ ^ /■ /■ /•y /•r / 

earth 


(b) 


segment above the ear 
the earth surface . 
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which the current and charge distributions can be assumed to be 
located along the wire axis, and the boundary conditions can be 
applied on the surface of the wire. 


At the surface of the conductor the sum of Incident electric 
field and scattered field is zero. Hence, 


E®(a) 


eJcs) 


(a) on wire surface. 


(3.2) 


where subscript 1 means that the component is in the direction of 
the wire axis. From (3.1) and (3.2) one can write 


E^(a) = j6iA^(s) + (^(s)) 


(3.3) 


where 


A^(a) 

<fr(s) 


I 

J axis 

J 

J axis 


I^(l’) cosd (a/l’)dl’ 
q(l’) G^(s/T’)dlV 





(3.4) 


(3.5) 


In (3.4) and (3.5) the Green's function and G* are given by 


G“ (a/T') 


u e 
*^0 


-Jk R 
•'os 


4nR 


(3.6) 


G* (s/1' ) = - 


-Jk R 
o a 


4nR £ 
a o 


(3.7) 


where R = I a - 1 ' 
a ' 


(3.8) 
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The current and charge density q(l’) are related by the 
equation o£ continuity, 

= TS JP- tii(i')) . (3.9) 

d is the angle between the directed length elements (axial 
directions) at the source and the observation points as shown in 
Fig. 3. 2. The source point 1' lies on the axis, and the field point 
s la on the surface of the conductor. 

Equation (3.3) is an integral equation for the potentials 
Aj(J)and^(I). Thld ia a mixed potential integral equation 
formulation in terms of current and charge distributions. In a 
mixed potential integral equation formulation a derivative of 
first order is present while in a single potential formulation a 
second order derivative la necessary. This gives the former 
numerical advantage over the latter. Since, charge and current are 
related by (3.9) the equation can be solved for current 
distribution in the wire antenna. The same solution procedure ia 
applicable whether the antenna ia in free space or it is above a 
lossy ground, the only difference being in the use of appropriate 
Green’s functions. In the next section the Green’s functions for 
the wire antenna in the presence of the dissipative half apace are 
discussed. 
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3.2.2. Wire antenna over lossy earth 

For a X directed horizontal current element over lossy earth, 

the X component o£ the magnetic vector potential in the air region 

is given by (2.24). The scalar potential is given by (2.61). 

Hence, the Green's functions G_(r/r’)and G.(r/^’) for the vector 

A 

and the scalar potentials, respectively, can be written as 


-jk R 
o 

(j e pj e 

- — 


-Jk R, 
o 1 


oo -)2y^h 


4nR. 


p P e 

T?- T — J (Xr)\dX 

j2n J o 


(3.10) 


-jk R 

G^(r/r’ ) = 


dO 


-Jk R, 2 -j2^ h 

■'ol p 

L_ + ® — J (Xr)XdX . 

4nR.4i; J2n I ,_2 .,_2 o^ ^ 

1 o 


0 


k.y +k J'- 
1 o o*^ 1 


(3.11) 


In (3.10) and (3.11) R and R^^ are the distances to the field point 
from the source and the image points respectively, and r is given 
by 


r = (X - x^)^ + <y - y^)^- (3.12) 

The points r and r' are respectively on the surface and on the 
axis of the wire. Further, it can be identified that r and r’ 
respectively stand for the points s and T’ in expressions (3.2) to 
(3.9). 
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From the Green’s functions given by (3.10) and (3.11), in 
principle any horizontal wire antenna above or on the ground can 
be analysed. However, numerical evaluation of these integrals is 
difficult and time consuming as the Bessel functions are highly 
oscillatory and converge very slowly. An accurate method for 
evaluating these Sommerf eld-type integrals is necessary as one has 
to compute these integrals repeatedly, for solving antenna current 
in frequency domain. For transient analysis the computation is to 
be repeated for many frequencies which makes it even more time 
consuming. Efforts to develop new, more efficient and more widely 
applicable methods for evaluation of these integrals are 
continuing [13,14,33,34,35]. 


Greenes functions with h ■ 0 


In the present investigation, our interest la about a loop 
antenna which la laid on the surface of the earth and hence, h the 
height of the antenna on the ground surface is zero. By putting 
h = 0 in Eqn. (3.10), the Green's function for magnetic vector 
potential can be written as 


G^(r/r’) 


h=0 



00 



Jo(Xr) 


XdX . 


(3.13) 


This Integral has been obtained in a closed form in Chapter 2, 
and la given by (2.38). Uniting g^(r/r’) for Green’s function for 
this special case (h = 0) one caJ^ write 
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e^(r/r ’ ) 


2nr(k^-k^) 
1 o 


a 


e 


-Jk,r1 

6 

T 


(3-14) 


From (3.11) the Green’s function for the scalar potential for h=0 
becomes 


G^(r/r’) 


2 

u 

< 

j2n 


00 


h=0 


1 


J^(Xr)XdX 

72 772 

‘"l^'o'^Vl 


(3.15) 


The Integral In (3.15) does not have a closed form expression. 
However, by putting k^ = 0, the Integral can be evaluated. Doing 
this amounts to obtaining fields in the air region by Laplace’s 
equation instead of Maxwell’s equation. This is permissible if the 
distance between the source and the observation point is very 
small in comparison to the free apace wavelength, or k^r << 1. 
This la a quasl-atatlc approach and will be valid for a practical 
sized loop antenna that may be used for underground mine 
coramuncation. In such cases it happens that |k^r| is quite large 
indicating that there may be considerable variation of current 
distribution in the antenna, but the condition k^r << 1 , is still 
satisfied. Writing k^(F/r ’ ) for the quasi-static Ck^^^O) 
counter-part of G^(r/F’) with h = 0 one gets 
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k^Cr/r') 


2 

j2n 


00 


J 


J CX.r)XdX 




2 

61 


2nk 


00 


I 


(Xr)dX 


1 0 


C3.16) 


2 

61 


2nk 




C3.17) 


3.3 POWER SERIES EXPANSION OF g.Cr.^»5 

A 


The Green's function 8^Cr/r’) given by (3.14) can be 


expressed in a power series form as follows. 


-Jk r ,2 ,.3 2 .43 

o . krjkr kr 

} 1 ., 

r “?“-^o~2l 3! 4! 


4 

jk r 
o 

51 


, kjr jkjr^ kjr^ 

'“F- = ? - - - 2 T — arr " -TT- 


si,5 4 
jk^r 

5 ! 


& 

av 


r -jk^r 
e e 


Jk^r I 


(k^-k^) JCkJ-k^jr 


1X1^ 1 ^ N ^ 

f (kj-k^)r 
8 


J(kj-k^)r^ 

Jo 


(3.18) 
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Multiplying (3.18) 
gets 


by 


2nr(k^-k^) 


and comparing with 


(3.14), 


one 


fi^(r/r') 





2j(kJ-k^) 

3(kJ-k^) 

1 o 


4 4 

(k, -k )r 
'• 1 o 

it ^1 2 , 2 . 

4(k, -k ) 
1 o 


15(k 



)r 


2 


'v 


■ f > 


J 


which can be expressed in the following form : 

OQ 

n=0 


(3.19) 


(3.20) 


where, 


^ = - 


o 

4n 


n 


2(-j) 
nl (n+2) 


(k 


n+2_j^n+Z 


) 


(k^ - 




(3.21) 


The expression (3.21) is not convenient for computational purpose 

because when k^=k , the numerator and the denominator both become 

zero. The factor (k.-k ) in the denominator can be removed as 

X o 

follows . 


n+2 n+2 

The term (k^^ ~ k^ ) appearing in the numerator of (3.21) can be 

factored by using the relation 


k” - k“ = (k. - k ) f 

1 o 1 o m— 1 


(3.22) 


where, £ 


m-1 


“o * •'o ♦ 


1 , m-2^, m~l 
k^ k +k 
1 o o 


(3.23) 


(m terms) 
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Removing the factor from the numerator and the denominator 

one gets 




n 


o 

Tn 


n 


2(-j) 
nl (n+2) 


n+1 


Cki + k^) 


(3.24) 


For the special case when k. 


one can write 


m-l 


= mk 


m-l 


k, =k 
1 o 


Hence, one gets 


n+l 


= Cn+2)k 


n+l 


k. =k 
1 o 


From (3.24) and (3.26) one can write 


(3.25) 


(3.26) 


n 


k, =k 
1 o 


_o 

4n 


2(-j) 


n (n+2)k 


n+l 


nl (n+2) 


2k 


And, hence 


Q 

4fi 


n I 


(3.27) 


fl^(r/r' ) 




o 

4n 


00 , . - n, n n-1 
(-J) k^ r 


I 


nl 


n=0 


(3.28) 


which can be identified as the power series expansion of the free 
space Green’s function. 
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A few of the coefficients given by (3.24) are as follows 


Q 

4n 

4n 


o 

4n 


j2 


/J 

Tn 


o 

4f7 


3Ck,+k ) 
1 o 


-1 


4(k +k ) 
1 o 


I k?+k k +k^ 1 

X 1 1 o o / 

I k^+kfk +k,k^ + 1 

\llolo oj 


- • T-g yi -T i. V- I kj+kfk + ... + k^ 1 

n 15(kj^+k^) \ 1 1 o o J 


4n • 72(k 


1 I 4 

^+k^) • t *"1 


k + . 
o 


+ k 




o 

4n ' 42 


— =r • i kf+kfk + ... + k^ I 

0(kj^+k^) X 1 1 ® ° J 


2880Ck 


1 / k^+k^k • 

k^+ko) I 1 1 o 


+ k 


) 


(3.29) 


o 

4n 


a26ftO(k 


rk-nry ■ { ■'X 


k + 
o 


. + k' 


} 


3.4 SOLUTION PROCEDURE 

For obtaining numerical solution to equation (3.3), the wire 
antenna is divided into N equal segments along the wire axis as 
showen in Fig. 3.3. The nth segment is indentified by its mid 
point n, and two end points n and n. The direction of the current 
is from n to n. Each segment Is sufficiently small, so that for 
any given segment the current can be treated as constant, with 
negligible error. 
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The total charges associated with the current Ij^Cn) in the 

nth segment of length Al^ are 4^ Ij^(n)/j<o. The positive charge is 

assumed to be distributed uniformly along the axis from the point 

n to a point at a distance of A1 along the wire axis. The centre 

n 

of this charge distribution is at point n and the linear apace 

which is being occupied by this distribution is identified by 

Likewise, for negative charge, n is the centre and the space is 

identified by A1-. Fig. 3.4 illustrates this model of a segment, 

n 

for current and charge distributions [32,39]. 


By approximating the derivative in (3.3) by a finite 
difference form and considering the mth segment, one can write 


^ ^ ^ VCm) - V(m) 

Ej^Cra) jwAj^Cm) + 

m 


(3.30) 


where, Aj^(m) is the vector potential at point m along the 
direction of the axis and is due to currents in all the segements; 
V(m) and V(m) are the potentials at points m and m due to the 

•f* 

charges associated with these currents. To obtain Aj^(in), V(m) and 
V(m), the Integrals (3.4) and (3.5) are to be evaluated along the 
entire wire axis. They are approximated as sum of integrals over N 
segements. Thus, one can write 


Aj^(m) 




g^(m/n)dl 



Ij(n) 


( 3 . 31 ) 
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(a) (b) 


Fig. 3*3 (a) A wire divided into N segments . 
(b) The nth segment. 



I 


Fig.3.4 Uniform current and charge distributions 
for the nth segment. 
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VCm) 


VCm) 


N 


!.{ q(n) J k^Cni/n)dl*+ qCn) J k^Cm/n)dl'J- 

n=l A1+ 

n 


C3.32) 


N 


qCn) J k^Cro/n)dl’+ q(n) J k^(m/n)dl ' j- , 

n=l A1+ 

n 


(3.33) 


where, cos & is the direction conaine between segements n and 
mn 

m, and q(n) and qCn) are the uniform charge densities with centres 

*4" “ "t* 

at n and n respectively. The charge density q(n) and q(n) can be 
expressed in terms of current I^Cn) as 


q(S) 


j<oAl + 


-I (n) 
qCn) = — 

joiAl- 

n 


(3.34) 


Equation (3.30) can be identified as an equation for N unknown 
currents. By allowing m to take all the values from in=l to m=N, a 
set of N equations can be obtained. These equations can be written 
in matrix form as 


[V] 


III [I] 


(3.35) 
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wh ere, f V] 


III 


eJcd-ai^ 

eJc2).A12 


tn 


I^Cl) 

IjCZ) 


eJcn) .A1^ 


IjCN) 


^11 ^12 ^IN 

^2 ^2N 


C3.36) 


-N1 


'N2 


'NN 


The equations defined by (3.35) can be viewed as equations for a 
N-port network. The matrices [V], [Z] and [I] stand for the 
network parameters, namely, voltage, impedance and current 
respectively. The element of the impedance matrix [2] can be 
written from (3.30) - (3.36) as 


Z 


mn 




where , 



mn 


jtAiCoa& A1 

mn ns 



g^(m/n)dl ’ 


k.(m/n)dl’ - f k.(m/n)dl’ 

* jcaAl- JaI- ^ 

+ ■' n n 

n 


— - — r k.(m/n)dl* + — - — j k (m/n)dl ' 

j6>Al+ * j6iAl- JaI- ^ 

^ n A1+ n n 




mn 


jcjAl 




(3.37) 

(3.38) 


(3.39) 
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Defining the integrals in (3.38) and (3.39) by 


^(ro.n) = 


1 

A1 


-■L 

n 


g^(in/n) dl ’ 


(3.40) 


?(in,n) 


— r 

A1 

-^Al 

n 


k.(in/n) dl’ 
<P 


(3.41) 


one can write 


= jwAl^.Al^'Km/n) + ^ [? (ro/n) - ?(m/n) - if(m/n) + ?(m/n)] 


(3.42) 


where A1 .A1 = A1 .A1 cos & 

n in n m mn 


(3.43) 


This result applies for self impedance (m = n) as well as for 
mutual impedance (m n). The integrals 4'(in,n) and f(m,n) are 
evaluated in the next section. 

3.5 EVALUATION OF 4' Cm»n3 


To evaluate ^(m,n), a local co-ordinate system with origin at 
n, and the y axis along the axis of the nth segment is constructed 
as shown in Fig. 3.5. From (3.20) and (3.40), iKm.n) can be 
written in the form 


A1 /2 
00 n 


^ I’’" 1 

n A— I J 

n=0 


-A1 /2 
n 


4'(m,n) 


dy’ 

m 


(3.44) 
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where# R 


m 


■/ p^+Cy-y ’ 


/ 2 ^ ,2 

ya +y’ 


TO n 


TO = n 


(3.45) 


By TOaking a change of variable 


t = y - y’ 

Al 


(3 . 46) 


and writing « 
rewritten as 


n 


, the integration given by (3.44) can be 


00 


*C».n) = J ^ b^I 


n 


(3.47) 


n=0 


y+a 


where# I 


= r fA^l 

n J [ J 


n-1 


dt 


(3.48) 


y-« 


The expression (3.47) can be written as sura of two series 

1 


‘Km.n) =2^ [S^ S^] 


where 


00 


I 


b I 
n n 


(3.49) 


(3.50) 


n=2i 


00 

I 

n=l+2i 


b I 
n n 


(3.51) 


with , i = 0 , 1 , 2 , 3 


m m M • f 


Expanding the series one can write 
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Fig.3.5 Geometry for evaluating iy^(m,n) 



Fig.3.6 (a) Circular loop antenna with excitation at g° 

(d) The reduced structure after being segmented. 
(Only eight segments shown ). 
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y+« 


J 


.2 . . 2 , 2 


= b I (t“ + p") "■ dt + b„ 

m 2 


y-« 


y+a 

J 

y-a 


.. 2 . 2.2 .. 
(t +p ) dt 
tn 


+ b, 


y+« 

J 

y-c* 


1 +— 

2 , 2 


dt + 


( 3 . 52 ) 


and Sj^ = ^3^3 ^ ^5^5 ^ 


y+a 


J 

y-a 


bjtbj I (t^ + p^) dt + bj 


y+a 

J 

y-a 


.. 2 ^ 2.2 .. ^ 

(t +Pjjj) dt + 


( 3 . 53 ) 


(a) Evaluation of S ; 

o 


The Integral ia given by 


y+a 

I 

y-a 


^^2, 2. 2 .. 
Ct +/»„) at 
in 


lOfi 


y+a+r p^ + (y+a)^ 
mi 

y~a+-/ p^+(y-«)^ 


( 3 . 54 ) 


The Other integrals of the series can be obtained by using the 
following recursive relation C36]. 
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I 


7 7 '^'^'7 

Ct^+ 2 dt 

m 


.^.2^ 2.^ 2 
tCt +p ) 


2(Y) + 1) 2(y) + 


C2T) + l)pf r ^2 

(t +p ) + C. 

1 


CJ p 

- f 

1) J 


C3.55) 


where y) = 0,1, 2, 3, ..., and is an arbitrary constant, 


(b) Evaluation S, 


The integral of the series can be written as 


2y )+ 1 


y+« 

1 

y-« 


(t^+p^)’^ dt, y) = 0,1,2 


(3.56) 


This integral can be evaluated by employing binomial expansion. 
One can write (for y? 0), 

^ ^ Oc Zt^Cn-i) ^ »cyt2("-2> + ... + p2„ 

m 1 m 2 m m 


(yj + 1) terms 


(3.57) 


Integrating (3.57), one gets 


y+ot 

J 

y-oi 


(t^+p^)^ dt 
m 


2yy-Hl y, 2^2(y,-l).l y, 4^2(yy-:t) + l 

, 1 «» . 2'^m . ,2y) 

2 y 7 + l ,a^-l) + l 2(r/-2) + l m 


y+or 


Jy-Ct 


(3.58) 


Following the above evaluation procedures for and , a 

few integrals (given by (3.48)) are as follows : 


-('-rSAL LiSRARV 

’ T-. KAiyplJR 


». No. 
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Io= log 


y+ot +y}r" 

m 

y-CA+-/Y” 

m 


I. = 2a 


y+a 


2 in 


vr^ yy- + -“ I 
2 ^ m 2 o 


1 nS) Zopf, 


1+4 1+4 o 2 

y+" Y 2 Y ^ + - 0^1 

"4" m “ m 4 <^10^2 


2/^>. 


1=4 YC5) + yC3) + Zap 

5 o 


(3-59) 


2 +Tr 2 +*^ c « 

y+g jr 2 _ y-g y 2 . 5 .2 
6 *in 6 m 


+ ^ P I/. 
6 4 


3p 

^ Y(7) * Y(5j' + P* YC3) + 2PP* 


where 


X = + (y+g)^ 

m in 


Y = + (y-g)^ 

m in 


(3.60) 


Y(n) = “ Cy“")”* n 1® integer. 

®Cin,n) is to be computed by using (3.29). (3.47). (3.59) and 

(3.60). 
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Evaluation of <;Cm,*n3 

Comparing (3.17) with (3.20), and (3.40) with (3.41), it can 

be seen that C(rofn) involves the integral I . One can write 

o 

2 

C(in,n) == 1 I. . (3.61) 

2nk2 

3.6 INPUT IMPEDANCE AND CURRENT DISTRIBUTION 

The current in the nth segment can be obtained by solving the 
equation (3.35). The solution can be written in the 
following matrix form. 

[1] - [Z]"^ CV] (3.62) 


It is assumed that the loop is excited by a unit voltage source at 
the gap which is located at ^ = 0*^ , as shown in Fig. 3.6. This 
excitation voltage la assumed to be applied to segment number 1. 
The voltage matrix can be written as 


1.0 



(3.63) 


rhe matrix [V] is a column matrix with N entries. The input 
impedance and the input admittance are given by 
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'in 


1 ( 1 ) 


in 


1 ( 1 ) 


(3.64) 


3.7 FIELD INSIDE THE EARTH 

The closed form expressions for potentials used to obtain 
current distribution in a loop antenna on the earth surface, 
cannot be used for computation of fields inside the earth. These 
expressions are valid only at the air/earth interface. To compute 
field inside the earth Sommerfeld type integrals representing the 
potentials in the earth region are to be evaluated numerically. Ue 
shall be computing magnetic field components only. The magnetic 
fields may be picked up using a small coil, preferably wound over 
a ferrite core, and then fed to the miner’s receiver. 


For a X directed current element, Fig. 2.1, the magnetic field 

components in terms of the magnetic vector potentials are given by 

(2.5a) ~ (2.5c). These components are denoted here as 

and <5H . to stress that these are elementary field components due 
zl 

to an isolated current element, and not the total field for a 
given antenna. These elementary field components can be obtained 
as follows : 


From (2.5a) and (2.27) one can write 
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<5H 


xl 




a 

dy 


CO 


_p S [■ 6xp{ jCr^z-?^^ 2 ^)}dX 

2n ^x I 2 ^ , 2 

*’ n k.y + k V. 

0 V o o*^ 1 


Sy 


px f’ (Xr)X^exp{ JC^^z-iKqZ^) )clX 


Zn 




,-1 2 , , 2 - 
rCk^j^o 


QO 


px r exp{j(y^2-y^z^))dX 


Zn 


J 


i 2 ^ . 2 

^ 1^0 + 


(3 


= K' 


where F^CXr) = ^ |J^(Xr)r“X - 2Jj^(Xr)/r 


) 


(3 


From (2.5b) one can write 


6H = 1- ^ (6A , ) - — J- (<5A , ) 

yl u k1 u Sx zl 


C3 


where the vector potentials are given by (2.26) and (2.27) 
(2.26) one gets 


1 ^ /-JIA s 

p az ^ ^xi^ 
'^o 


a 

Sz 


p 

2 n 


00 


P 

j 2 n 


J 


J^(Xr)Xexp{ j(yj^ 2 -^'^ 2 ^) )dX 




00 


Ic 


jQ(Xr)^'^Xexp{ j(^'j^z-j'^z^) } dX 
_ __ _ 


(3 


.65a) 


. 65b) 


. 66a) 

From 


.66b) 
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From C2.27) one gets 


<3Cl 




li 


A 

2n 


Sx 


m 


Cr 1 (Xr )\ exp { j Cy 2 2-y^z:^ ) } dX 


1,2 2 


(3.66c) 


2n 


J 


^ Cyj-y^)J^CXr)X^Fy(Xr)exp{ 


1.2 2 

k.,y +k y. 

1 o o 1 


where F (Xr) 

y* 


J^(kr) 


J{1^ 


(Xr)X 


2Jj^(Xr) 

r 



(3 . 66d) 


C3.66e) 


From (2.5c) and (2.26) one can write 

= {h 



oo 

(Xr)X exp{ jCj'j^z-^'^z^)} dX' 

1 

T^Jo" 


„ fy 

JjCXr)X 

2 

03Cp{ j(yjZ-^^z^))dX 

iZnv 

1 

0 

y .< +y 

1 o 


These integrals for the elementary field components are oscilla- 
tory and converge very slowly. Uhen the observation point is 
inside the earth the factor exp { j (y ^^z-y^z^) } enhances convergence 
of these integrals. The convergence is quite fast when the 
observation point is located at a considerable depth, which is 
usually true for underground mines. To compute the field 
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components for the whole of the antenna on the earth surface, 
these elementary field components are to be integrated for all the 
current elements along the wire axis. 

3. 8 NUMERICAL RESULTS 

Before doing actual computation it is necessary to ensure 

that the quasistatic approach i.e. , k b << 1, is valid for the 

o 

loop antenna being analysed. The purpose of these antennas, is to 

transmit voice signals to underground miners. For a voice grade 

communication the highest frequency is 3000 Hz. In the present 

study the maximum size of the loop is considered to be 600m. Uith 

these values we have the maximum value of 2k b = 0.038. This value 

o 

of 2k^b << 1, and hence the quasi-static formulation for the 
scalar potential can be used. 

In figures 3.7 and 3.8 the variation of the real, R (I(<^)), 

© 

and imaginary, I parts of current are shown. For a 

m 

small loop, I (IC<ii3) remains almost constant for the whole loop; 
m 

for example for the 50m loop I^(IC0°)) = 165.77 mA and 

I (1(180*^)) = 165.63mA. This is not true for R^(ICi^l) which can be 
m ® 

seen from the fact that R (I(0'^))=8.67mA and R (I(180°)) = -3.84mA. 

6 “ 

This variation of R (I(<ifr)) in a small loop on the earth surface or 

e 

in a homogeneous lossy medium makes it different from a small 

loop in free space where R (I(<J^)) remains almost constant. 

© 

However, since for a small loop I (I(<^)) >> R (I(<;6)), the loop 
current can be approximated by 1^(1 (<&)), which is the usual 
criterion for identifying a small loop. As the loop diameter is 
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increased R^C I (<?!>)) becomes comparable to and hence the 

magnitude and phase of the current vary considerably with 4> * which 

can be seen in Figs. 3. 9 and 3.10. Uhen b is considerably large the 

current drops very fast near the feed region of the loop. For 

example, for 2b = 100m, the ratio j IC180°) j / | 1(0°) j is 0.97; 

whereas for 2b = 200ra and 300m, these values are 0.70 and 0.39 

respectively. A large value of this ratio indicates that the far 

end section of the loop (diamet erically opposite to the feed 

region) has practically no effect on the feed region. This effect 

can be seen in Fig. 3.11, where the variation of input conductance 

G and susceptance S, with the diameter 2b is shown. Both G and S 

remain almost fixed when b is large. Uhen b is small G varies 

linearly with o’, and G is much smaller than jSj. When b is 

increased, the effect of o’ becomes more and more prominent, which 

finally leads to the saturation for the input admittance values. 

Similar trends are also observed when a loop is fully embedded in 

a conducting medium [9]. In Tables 3.1 to 3.9 the input 

admittance Y along with 1(0°) and 1(180°) for different values of 
n 

parameters have been tabulated. From these tables it can be seen 
that when |Sj is at least seven times greater than G, the ratio 
R = 1 1 (180°) / 1 1 (0°) I is greater than 0.95. For practical purpose, 
a loop antenna having greater than this value (0.95) for can be 
considered as a small loop. 

Figures 3.12 and 3.13 show the magnitude of magnetic field 
components inside the earth due to surface current distributions 
in the loops. The field strength has been obtained by normalizing 
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the feed point current 1(0°) to unity. Fig. 3.12 is for the case 
when the current distribution is almost uniform, while Fig. 3. 13 is 
for when their is a considerable variation of current distribution 
with These current distributions are shown in Figs. 3. 7 to 3.10. 
The fields depend both on radial distance from the loop axis and 
angle When the current distribution is almost uniform the H 

z 

component practically does not vary with while component 

varies to a small extent. The situation is not the same when there 

is a considerable variation of current distribution. In Fig. 3.13 

it can be seen that when <^ = 0° the radial component is much 

larger than H or H values at any other angles. This is due to 
I? z 

the fact that current is mainly concentrated near the feed region 
of the antenna. The peak occurs approximately straight below the 
feed point. 

3.8.1 DISCUSSION ON THE NUMERICAL METHOD 

The formulation of the matrix equation (3.35) is based on the 
assumption of an uniform current for a given segment; and using a 
single matching point for the scattered field for each segment. 
This procedure is equivalent to using rectangular basis and delta 
test functions. Use of these simple functions requires 
sufficiently large number of segments to approximate the current 
distribution so that the current in two neighbouring segments does 
not change appreciably. However, the major advantage gained by 
using these functions in conjunction with the closed form 
expressions for the potentials, la that the solution procedure 
does not require any numerical integration for obtaining current 
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Fig.3.7 Variation of the real part (Rc(I) ) of the loop 
current with 0 . 




I m (I) with 0 . 



Phase angle ( radians) 





Phase angle (radians) 





Susceptanca , B (mir) Conductance ,G (mv) 
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Fig. 3.12 Magnitude of magnetic field components Hz(^) and 
Hr(0)inside the earth due to current distribution in the 
surface loop. 


Hr (I?) or 
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TABLE 3.2 


Input (xcJ-Tfi (Itance for a circular loop antenna on the earth surface 


f = 

1500 Hz, o- = 0.001 

mho/m, 

11 

a = 0.0025 m 


Loop dia 

in 

I (o‘ 

^) * 

I (180 

°) * 

(at) 

(mini mho) 





50 

4.86- j3<54.07 

364.1 ^ 

/-I. 56 

364.1 / 

-1.58 

100 

8.89-jl69.48 

169.7 ^ 

/-I. 52 

169.4 / 

-1.59 

150 

12.74-jl08.80 

109.4 ^ 

f_-l . 45 

108.5 /_ 

-1.62 

200 

16.46-j79.85 

81.5 ^ 

/-I. 37 

79.0 ^ 

-1.66 

250 

20.03-j63,29 

66.4 ^ 

/-I. 26 

61.8 ^ 

-1.72 

300 

23 . 46- j52 .92 

57.9 ^ 

/-I. 15 

50.4 ^ 

-1.78 

350 

26.69-J46.14 

53.3 ^ 

/-I. 05 

42.1 ^ 

-1,85 

400 

29 . 69- j41 . 68 

51.2 ^ 

/-0.95 

36.0 1 

-1.94 

450 

32 . 40-J38.81 

50 . 6 ^ 

/-0.86 

31.1 ^ 

-2.03 

500 

34.76-j37.09 

50.8 ^ 

/'-0.82 

27.1 ^ 

-2.13 


* Magnitude in inA and phase angle in radians. 
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TABLE 3. a 


Input O-d'Ti'ntJanc 6 for a circular loop antenna on th.6 earth surface 


f = 

1500 Hz , O' = 0 . 001 

mho/in, £ 

= 15, 
r 

a = 0.0025 m 


Loop dia 

y, 

in 

1 (0‘ 

* 

I (ISO 

A 

(m) 

(fflilli mho) 





50 

4.86- j364.07 

364.1 ^ 

/-I. 5 6 

364.1 / 

-1.58 

100 

8. 89- jl69 . 48 

169.7 ^ 

/-I. 52 

169.4 / 

-1.59 

150 

12 . 74- jl08.80 

109.4 ^ 

/-I. 45 

108.5 / 

-1.62 

200 

16.46-j79.85 

81.5 ^ 

^-1.37 

79.0 ^ 

-1.66 

250 

20.03-j63.29 

66.4 ^ 

/-I. 26 

61.8 / 

-1.72 

300 

23.46-j52.92 

57.9 ^ 

/-I. 15 

50.4 / 

-1.78 

350 

26.69-j46.14 

53.3 ^ 

/-I. 05 

42.1 / 

-1.85 

400 

29 . 69- j41 . 68 

51.2 ^ 

/-0.95 

36.0 / 

-1.94 

450 

32.40-J38.81 

50.6 ^ 

/'-0.S6 

31.1 / 

-2.03 

500 

34.76-j37.09 

50.8 ^ 

/-0.82 

27.1 / 

-2.13 


* Magnitude in mA and phase angle in radians. 
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TABLE 3.3 


Input o.d'nnilla^iC'S for a circular loop antenna on the earth surface 


f = 

1500 Hz, O’ = 0.002 

mho/tn. 

£ =15, 

r ' 

a = 0.001 m 

Loop dia 

in 

I (0 

* 

I (180°) * 

(m) 

(milli mho) 




50 

8.56- j332 .84 

332.9 ^ 

-1 . 54 

332 .7 / -1.58 

100 

16.14-jl55.48 

156.3 1 

-1 . 46 

155.1 ^ -1.62 

150 

22.95-jl01.20 

103.8 1 

-1 . 35 

99.9 / -1.68 

200 

29.88- j75. 28 

81.0 

1.19 

72.3 ^ -1.76 

250 

36 . 10- j61 .60 

71.4 1 

-1 . 04 

56.1 ^ -1.86 

300 

41.68-j54.05 

68.3 1 

-0.91 

45.1 1 -1.98 

350 

46 . 41- j50 . 18 

68.4 1 

-0.82 

37.0 ^ -2.12 

400 

50 . 12- j4S . 61 

69.8 1 

-0.77 

30.6 / -2.28 

450 

52.78-j48.42 

71.6 1 

-0.74 

25.3 / -2.43 

500 

54.43- j48. 93 

73.2 1 

-0.73 

20.8 -2.60 


* Maenltude in mA and phase angle in radians. 
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table 3.4 


Input adwMiance £ot a circular loop antenna on the earth surface 


£ = 

1500 Has, & = 0.002 

mho/m, 

% = 15, 

a = 0.0025 

m 

Loop dia 

Y 

nn 

I 1 

(0°) * 

I 

(180°) * 

(ffl) 

(ntilll mho) 






50 

9 .64-j364.11 

364.2 

^-1.54 

364 . 1 

z 

-1.58 

100 

17.74“J169.75 

170.7 

/-1.46 

169.3 

z 

-1 .61 

ISO 

25.40-J109.67 

112.6 

^-1.34 

108.2 

z 

-1 . 67 

200 

32.64-j81.76 

88.0 

/-I. 19 

78.6 

z 

-1.75 

250 

39.35“j66.80 

77.5 

^-1.03 

60.8 

z 

-1.85 

300 

45.36- j58. 56 

74.1 

^-0.91 

48.9 

z 

-1.98 

350 

50.45-J54.31 

74.1 

^-o.a2 

40.0 

z 

-2.12 

400 

54.44~J52.56 

75.7 

/-0.76 

33.0 

z 

-2.71 

450 

57.29-J52.31 

77.6 

/-0.74 

27.3 

z 

-2 . 43 

500 

59.06-J52.83 

79.2 

^-0.73 

22.4 

z 

-2.59 


* Magnitude in roA and phase angle in radians. 
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TABLE 3.S 

Input fot a circular loop antanna on the earth surface 


£ = 

3000 Hz, a = 0.001 

Miho/in, £ =15. 

r ' 

a = 0.001 Ml 

Loop dla 

In 

I (O'^) * 

I (180°) * 

(Ml) 

(willi Miho) 



50 

4.34- jl65.70 

165.8 /-I. 54 

165.7 / -1.58 

100 

8.07-J77.72 

78.1 ^-1.47 

77.5 ^ -1.62 

150 

ll.60-j50.36 

51.7 /-I. 34 

49.7 ^ -1.68 

200 

14.94-j37.61 

40.5 ^-1.19 

36.2 / -1.76 

250 

18.06- j30. 77 

35.7 ^-1.04 

28.1 ^ -1.86 

300 

20.85- j26. 99 

34.1 ^-0.91 

22.5 ^ -1.98 

350 

23.22-j25.05 

34.2 ^-0.82 

18.5 ^ -2.12 

400 

25.08“ j24. 27 

34.9 /'-0.77 

15.3 / -2.28 

450 

26.41~j24.17 

35.8^-0.74 

12.7 /-2.44 

500 

27 .24“ j24.43 

36.6 ^-0.73 

10.4 / -2.60 


Magnitude in tnA and phase angle in radians. 
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TABLE a.e 


Input cud rdi tilance £or a circular loop antenna on the earth surface 


f = 

3000 Hz, O' a 0.001 

mho/ffl, = 15, 

a = 0.0025 m 

Loop dia 

(m) 

Y 

in 

(mini mho) 

I (0°) * 

I (ISO'^) * 

50 

4.82“jl82.05 

182.1 ^-1.54 

182.0 ^ -1.58 

100 

8.87-J84.86 

85.3 ^-1.47 

84.7 ^ -1.62 

150 

12 . 70” j54.81 

56.3 ^”1.34 

54.1 / -1.68 

200 

16 .32” j40.85 

44.0 ^-1.19 

39.3 / -1.76 

250 

19.69-J33.36 

38.7 ^-1.04 

30.4 / -1.86 

300 

22.69-j29.24 

37.0 /-0.91 

24.4 / -1.98 

350 

25.24-j27.11 

37.0 1 - 0.62 

20.0 / -2.12 

400 

27 . 25” j26.24 

37.8 /-0.77 

16.5 / -2.27 

450 

28.67-J26.12 

38.3 ^”0.74 

13.7 / -2.43 

500 

29.56-J26.38 

39.6 ^-0.72 

11.2 ^ -2.59 


* Maghltude in inA and phase angle in radians. 
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TABLE 3.7 


Input adiDilianee for a circular loop antenna on the earth surface 


f = 

= 3000 Hz, a ~ 0.002 

mho/m, 

£s ^ - 15 , 

a = 0.001 m 

Loop dla 

Y 

^In 

I (0 

°) * 

I (180°) * 


(mini mho) 




50 

8. 67- jl65. 77 

166.0 1 

-1.52 

165.7 ^ -1.59 

100 

16 .09-J78.21 

79.8 1 

-1.37 

77.4 ^ -1.66 

150 

22.96-j51.87 

56.8 ^ 

-1.15 

49.4 ^ -1.78 

200 

29 .08“ j40.87 

50.2 ^ 

-0.95 

35.3 ^ -1.94 

250 

34.07-J36.39 

49.9 i 

-0.82 

26.6 ^ -2.13 

300 

37 .58- j35.20 

51.5 1 

-0.75 

20.3 ^ -2.35 

350 

39.59- j35. 53 

53.2 ^ 

-0.73 

15.4 / -2.58 

400 

40.37- j36. 30 

54.3 1 

-0.73 

11.5 ^ -2.81 

450 

40.41-j36.95 

54.8 1 

-0.74 

8.1 ^ -3.03 

500 

40.11“j37.32 

54.8 1 

-0.75 

5.3 ^ -3.04 
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TABLE 3.8 


Input c.chT> iLtance for a circular loop antenna on the earth surface 


f = 

3000 Hz, a = 0.002 

mho/m, £ =15, 

r ' 

a = 0.0025 

m 

Loop dia 

Cm) 

in 

(mini mho) 

I (0°) * 

I (180°) * 

50 

9 . 61- jl82. 12 

182.4 /-I. 52 

182.0 ^ 

-1 . 59 

100 

17.69- j85. 40 

87.2 /-I. 37 

SI. 5 1 

-1 . 66 

150 

25.14-J56.46 

61.8 /-I. 15 

53.7 / 

-1 . 78 

200 

31.76-j44.40 

54.6 /-0.95 

38.4 / 

-1 . 94 

250 

37. 13- j39.47 

54.2 ^-0.82 

28.8 / 

-2 . 13 

300 

40.92-J38.13 

55.9 /-0.75 

22.0 ^ 

-2.35 

350 

43.07-J38.45 

57.7 ^-0.73 

16.7 l_ 

-2.57 

400 

43 . 91-J39.25 

58.9 /-0.73 

12.3 / 

-2.80 

450 

43.95“j39,94 

59.4 ^-0.74 

8.7 / 

-3.02 

500 

43.62-J40.34 

59.4 ^-0.75 

5.6 ^ 

-3.06 


Magnitude In mA and phase angle in radians. 
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table 3.9 


Input for a circular loop antanna on the earth eurface 


f = 

3000 Hz, O' = 0.01 

mho/m, = 15.0, 

a = 0.0778 

m 

Loop dia 

Y, 

in 

I (0*^) * 

I (180°) * 

(Ht) 

(mho) 




10 

0.025~jl.989 

1.989 / -1.56 

1.990 ^ 

-1.57 

20 

0.040-J0.857 

0.858 ^ -1.52 

0.856 ^ 

-1 . 59 

30 

0.054- jO. 529 

0.531 ^ -1.47 

0.527 £ 

-1 . 61 

40 

0.068-J0.378 

0.384 ^ -1.39 

0.375 /_ 

-1 . 64 

50 

0.0S0-J0.293 

0.304 / -1.30 

0.288 ^ 

-1 . 68 

75 

0.110-J0.192 

0.221 ^ -1.05 

0.177 ^ 

-1.81 

100 

0.134-j0.153 

0.203 / -0.85 

0.122 ^ 

-2.01 

125 

0.152-j0.139 

0.206 1 -0.74 

0.088 / 

-2 . 23 


* Magnitude in ampere and phase angle in radians. 
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distribution in the loop. All integrations have been expressed in 
closed forms. 

For numerical computation, the circular loop has been 

approximated by 36 equal and linear segments. To evaluate i'Cm.n) , 

the first eight terms of the series (3.47) have been used. Xo 

check correctness of the numerical method, circular loop antennas 

in free space have been analysed using similar series expressions. 

with eight term approximation, and these numerical results have 

been found to be in close agreement with the published results [9] 

f or k b SI 1.3, 
o 

5 

For a very small loop which almost carries an uniform 

current, the computation of input impedance by this numerical 

method is difficult to accomplish. A small loop (electrically) can 

be identified by a small value of <o. When « is very small the 
A 

Impedance 2^^ due to the magnetic vector potential is also very 

A A \ 

small, as ot <*i, but on the other hand, . Since, for 

small values of w, >> the impedance matrix [Z] = 

[Zf^+Z^^] becomes due to the limited width of memory in a 

inn mn 

computer. This matrix is singular, since an uniform current 
distribution does not depend on depends on To 

have a uniform current the effect of charges are to be cancelled 
mutually and it requires subtracting large and almost equal 
numb.ro ropr.oontlna Practically, thlo lo difficult to 

accomplioh in a ccput.r. In the pr.o.nt InyeotlBatlon th. 
comput.r proBrann. hao bean v-rltt.n in FORTRAN with doubla 
pr.clalon arlthnatic calculation.. It has bean found that a loop 
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of ten metre diameter for frequency 3000 Hz and o- = 0.002 mho/m 
can be analyzed. This corresponds to (k^^jb i 0.03. A loop with 
|kj|b < 0.02, does not lead to convergence. In Chapter 4 a closed 
form expression will be developed for self-impedance of a circular 
loop antenna on the earth surface. 


3.9 THE LOOP ON THE EARTH SURFACE AS A HOMOGENEOUS MEDIUM PROBLEM 


The loop antenna on the earth surface can be viewed as a 
homogeneous medium problem by using an effective propagation 
constant. The following effective propagation constant has been 
suggested by Popovic [37]: 





(3.68) 


It has been found that numerical results obtained by using k^.,, 

et t 

are marginally different from the results shown in Figures 3.7 to 
3.11, which have been obtained via Sommerfeld type integral 
representations. To understand this, a comparison of the 
potentials obtained from k^^^^ with the potentials used in this 
chapter will be useful. 


The vector and the scalar potentials due to a current element 
of length dl, at a point at distance r from the source can be 

written in terms of as 




e 


- jrk 


ef f 


(3.69a) 


4nr 
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qdle 


Jrk 


ef f 


4n r<o 


9ff 


where , 


u- 

e£f 

ef f 2 
<<o ju 


(3 . 69b) 


(3 . 69c) 


The expreaalons (3.69 a) and (3. 69b) can be written in the 
following aeriea form : 




fj Idl 


TfT 


(3 . 70a) 


eft 


(3.70b) 


where 


f(r,k^,kj^) » i 


(kjtkj)l/^ r(kjtkj) r^kj.k^)^/^ 

j + j_ 




12-/2 


3., 2^, 2.2 

r (k.+k ) 

— n + 


(3.70c) 


Expreaaiona (3.69a) and (3.69b) can be compared with (3.19) and 
(3,17). Since, k <<lkj, we can approximate these aeries 

O i 

expreaalona by putting k^=0. Uith this approximation one gets 


*^eff 




k! 


^ ff 


(3.71) 
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and the expreaeion (3.19) becowea 


fi^Cr/r ' ) 


k =0 
o 


«= !!£ [i 

4 n 't 


j 3"^! 


1 2 ,3 2 4 3 

k^r kjT 

T~ ■*■ j —TT — + — cn — +. 

4 15 72 


From the above one can aee that 

(1) the flrat term of the aeries (3.70b) is identical 
scalar potential given by (3.17) and 

(ii) the series representations (3.70a) and (3.72) match 
at least upto first few terms. 


.] . 

(3.72) 

with the 

closely, 



CHAPTER 4- 


SELF IMPEDANCE OF A SMALL LOOP ANTENNA ON TI-CI EARTH SURFACE 

4, A INTRODUCTION 

Ifi this chapter wa coiiaidar a aotall loop antenna on the earth 
aurlacje atid aaaiiflte that the loop la carrying an iinl£onii current. 
The aelf latpedance (input Itnpedance) o£ a aoiall loop In free apace 
can be computed by adding the radiation reaiatance with the aelf 
inductance (reactance) of the loop. The Inductance can be computed 
uaing Neumann' a formula. The radiation reaiatance la derived by 
integrating the radiated power over a large aphere. Uhen the loop 
la on the earth aurface, there will be dlaalpation In the body of 
the earth, aa well aa radiation to the free apace. The equivalent 
reaiatance which repreaenta theae power outflowa from the antenna 
can not be derived by taking the aame approach uaed for radiation 
reaiatance in free apace. The Neumann’a formula in ita baaic form 
la expreaaed in terma of a d.c. magnetic vector potential. Ue 
ahall replace thla d.c. potential by a time varying potential for 
computing the aelf impedance for a small loop antenna on the earth 

aurface. 

4, a NEUMANN* S FORMULA 

Fig. 4.1 deacribea a circular loop made up of a conducting 


wire. 

The mean 

diameter 

of the loop is 2b and the 

radius 

of 

the 

wire 

is .a. The 

self inductance of the loop is 

composed 

of 

two 

parts 

! One part 

is due 

to the magnetic flux 

enclosed 

by 

the 
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Fig. 4.1 A circular conducting loop. 
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contour and the other part i. due to the flux inside the 
conductor- Let is the inductance due to the flux enclosed by 
C^. This inductance can be obtained by using the Neumann’s formula 

4^ ^0-^1 ■ C4.1) 

This expression assumes that the current distribution inside the 
conductor can be replaced by a filamentary current defined by 
contour C to obtain the external flux enclosed by the contour C- . 
Our aim la to Incorporate necessary modifications to (4.1) so that 
it can be used for the computation of self impedance (resistance 
and reactance) of small antennas on the earth surface. 


4..2I.1 Modifleid N«eutnanni'’K Formula 


The right hand side of (4.1) can be identified as the 
Integral of the magnetic vector potential along the contour 
due to an unit d.c. current flowing along the contour C^. If the 
current Is time varying with a small angular frequency ei, the 
induced e.m.f. associated with is given by 




LI 


C. C, 


4 ^ 


dl -dl- 
o 1 


(4.2) 


Let us assume that the wire is a perfect conductor. Then. under 
time varying condition there will not be any flux inside 
conductor. The Inductance L^ is now equal to the self inductance 



of the loop. The applied voltage which must act opposite to E 
given by 
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Expression (4.3) ignores the phase change experienced by the 

magnetic vector potential under time varying condition. To take 

this Into account the d.c. magnetic vector potential is to be 

replaced by the time varying magnetic vector potential. For free 

fj - jk R 

space we can replace by p^e ^ /4nR, and for a small loop on 

the earth surface the replacement is to be done by the magnetic 

vector potential given by (2.38). Self impedance Z of a loop la 
equal to the applied voltage when the loop carries an unit 

current. Hence Z for a constant current carrying loop on the 

a 

earth surface can be written as 




2n(kj-k;) •'c^ 


o 


•jk^R 


jkj^R 


R 


dl .dl. 
o 1 


(4.4) 


Expression (4.4) can be written in a series form by using (3.19) 
as 


J 

* -tr^ le 

O 1 




dl .dl. 
o 1 


(4-5a) 
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whore 

« 1 


• 



B„ » -■ 

j2 

■'i 

- 

0 


2 

3 

1 1 

- k2 

oj 


r ^ 



(4.5b) 


Rotaibing only the first four terms of the series, one can write 


& 




I^(R) + B2I2CR) + B3I3CR) + B^I^(R) 


(4.6a) 


where 


Ij^(R) » 

/ <e) 

^0 '“i 

(4.6b) 

igCR) “ 

j j 

0 1 

(4.6c) 

I3(R) « 

j r RdT^.dXi 
<=1 

(4.60) 

I^CR) * 

I S «'3T„.drj 
<^0 <^1 

(4.6e) 
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4. a. 2 Evaluation of the integrals 

First we consider the Integral I^(R) The 
3X^ And are given by 

1 ^ 0 1 " , |dTj = (b-a) dd^ - b de^ 

The angle between dX^ and dTj^ la and hence 


3T .XT. » b^cosC© -d ) dd de, 
o 1 ' 1 o' o : 


The 


distance R between XT^ and dl^ la given by 


r 2 - b"* t (b-a)2 - 2bCb-a) coa(^^-^^) 


From (4.6b) and (4.7a) - (4.7c), one can write 

Zn 2n 


I,(E) = I I 


COB ) dd. d& 

1 o 1 o 




0 "0 |b^+(b-a)^-2b (b“a)cos(e^-e^)j 


Ue may Integrate over first (for a fixed value of 


(6»j^ ~ *“ ^ , oae gets 


2n Zn 


Ij_(R) » b ^ J I 


cos6f dd dd 


0 0 Jb^+(b-a)^ - 2 b(b-a)cosd j 


1/2 


magnitude of 

(4.7a) 

(4.7b) 

(4.7c) 

__ . 

(4.8a) 

d^') . Defining 


(4.Sb) 



This integral may be evaluated In any order. Integrating 
one gets 


2n 

« Znb^ 


where 

[ 2 9 n I / 2 

b + (b - a) - 2b(b - a) cosS 

“ la^ + 4b (b - a) ain^(d/2)j 


I 


Coa6> do 


The expreaaion (4. Sc) can be evaluated in terms of 
integrals [38]. The final result is 


Ij^(R) * 4nb 




2 


where « 4b <b - a)/(2b - a)^ and K and E are elliptic 

lit 

given by 


K » 


I 


n/2 


dx 


0 (1 - kl sin^x)^^^ 

49 


w/2 


E » 


I 


(1 


. 2,2 . 1/2 . 

k:sln x) ax 
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over & 

o 

(4.8c) 


(4.8d) 

elliptic 

(4.9a) 

int egrals 

(4.9b) 

(4.9c) 



87 


For a << b, the result o£ (4.9a) reduces to 


I^(R) = 4fib |ln(8b/a) - 2| 


(4.10) 


By golne through altnilar exercises which were performed to 

arrive at (4-8c) from (4.8a), the integrals I.,(R), I.,(R) and I. 5 (R) 

X 3 

can be written as follows. 


I_(R,) s 2n:b‘' 


2n 

J 


COB& d9 


m 0 


(4.11) 


2n 

I^CR) » 2nb^ J 


t1/2 


a^+4b (b - a) aln^(e/2) 


coa& dd 


2n 


4 

I 


(4.12) 


2n 


I^(R) « 2rxb^ J Rj cos© dd 


(4.13) 


Since, a«b, R^^ in (4.12) and (4.13) can be approximated as 


Rj^ tn 2b «ln(0/2) 


(4.14) 
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Substituting frow (4.14) (4 12 ) 

2n 


one can write 


I 3 CR) “ 4nb^ J coa 6 .aln(e/ 2 ) 


d<9 


= ~ 


Similarly one gets 


In 


I^CR) « 8 nb^ J cos©.ain^(e/ 2 ) 


dd 


0 
2n 

» 4nb^ j (1 - cosd) 

■^0 


“ -4n^b^ 


Substituting I^(R) where (1 . 1 . 4 ). in (4.6a) and using 
one can write 


■ ^.1 " ^,2 ♦ 2,3 * 2,4 


where 


^sl " ■{ln(8b/a )-2 


2 « « 0 
a2 


{■ 


) 


. . kj + kjk + k.k^ + k^ 

« „ t.3 / 1 1 o 1 o o. 

» JC-t) C iT-i-c ) 


‘‘1* “ o 


^a4 ** n6>^^b ( 


. k 5 + kjk + kjk^ + k k^ + k^ 
4.1 lo lo lo o 


*^0 


) 


(4.15) 


(4.16) 
(4.5b) , 

(4.17a) 

(4.17b) 

(4.17c) 

(4.17d) 


(4.17el 
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When the loop la In a homoeeneoua medium (loaay or free apace) 
with “ k, one obtalna 


'a3 


1 o 


J<§) ^ k^ 


(4 . 18a) 


2 


a4 


k. 


=k “k 
o 


C-?)" 


(4. 18b) 


With thla, the expreaalon for 2 becomea Identical to the 

a 

expreaalon filven In f9] for a amall loop in a homogeneous mediuoi. 

It can be seen that for free apace, Z . and Z „ are reactive. 

al a3 

vheneaa 2 ^ la resistive. The inductance associated with 2 . can 

8 4 S J. 

be identified as L the inductance of a loop due to the external 
magnetic flux. The term 2^^ is the radiation resistance for a 
small loop in free space. 


For a loop antenna on the earth surface and operating at low 

frequencies, |k. j >> k^. This allows one to put k^=0 in (4.17d) 

2 2 

and (4,17«) without any significant error. Since k^ =£ -Jib'll 
can write 


■ 40ftk^|kj| b^ 


(4 .19a) 
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Since a* (4.19a) can also be wrltt 


en as 


<7 ,1 2 2 . ^3 

^S3 I 


(4.19b) 


Thlfl tenm Is resistive. The contribution from Z is both 

Q 4_ •-r w fc. ♦* 

resistive and reactive, but is the dominant term. Neglecting 
the contribution £rom the self impedance can be written as 

^s “ ”^1 (4.20a) 

where 


Rj^ « 40nrk^|kjjb^ 


(4 . 20b) 


Xj^ “ 120nk 


<.•> (‘"<-1^) - 


(4.20c) 


It tan b. .e.n that for a omall loop 1!^ « b . Since, « 

X, , conductance , and ausceptance B , can be written as 

La a 

R, 


Q M 

a 


h ^ h 



and 


* 


-JX. 

4 * 


ftt " J 
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4.4 Dis:cussa:oN of results 

Comparison of results obtained for small loops (where |B^j is 
atleaat ten times greater than G^) by using expression (4.20) and 
by the method of momenta discussed in Chapter 3 leads to the 
following observations : 

(i) Susceptanee values are in close agreement with each other. 

(11) The conductance obtained using (4.20) is about an order of 
magnitude smaller than the values obtained from the method of 
momenta. For the particular cases considered in Chapter 3 
this difference has been found to be by a factor of 25 to 55. 

(ili)The conductance given by (4.20) does not depend on the radius 
of the wire, but the values obtained by method of moments do 
depend on It. 

The observations (1) to (ill) also hold good for loop 

antennas embedded in a hemogeneoua dissipative medium [9J. It is 

useful to determine using a simple calculation procedure how large 

a physical loop can be treated as an electrically small loop. Ue 

G 

consider a small loop as "small” when | ^ 0.15. Uhen this 

ratio increases beyond this value the input admittance reaches 
saturation very rapidly. Since, can be computed accurately 

using expression (4.20), it la desirable that the conductance 
(or resistance R^) la alao computed accurately using a simple 
formula. It has been observed that if the input resistance value 
given by (4.19a) or (4.19b) is multiplied by factor which la 

given by 



92 


where 


F 


c 


X 


2n 

3 


O 


a = 21n(-^) 
61 


(4.21) 


the new value agreea closely with that obtained by the method of 

momenta within 12 percent. Using the factor F one can write 

c 

input resistance for a small loop antenna on the earth surface 

and operating at low frequency (so that at -Jjkjj) as 


n 2n 2 2 ,3^ 

R , ® ^ o'b rt 

in 9 o 



(4.22) 


It may be pointed out that for a small loop (k^b < 0.2) in 
free space the radiation resistance given by (4.18b) can be used 
as the input resistance for a loop antenna. The disagreement 
observed for a loop on the earth surface is due to the fact that 
the real part of the current for a loop antenna on the earth 
surface varies to a large extent, while the imaginary part remains 
almost constant. Therefore, the assumption that the current is 
constant for a small loop on the earth surface, although can be 
made for the reactive part, the same does not hold good for the 
real part. 



CHAPTER 5 


PERFORMANCE OF COHERENT PSK RECEIVERS IN MINE ENVIRONMENT 

5.1 INTRODUCTION 

Th6 p ( otT'itt&iic 6 of A ELF/VLF cecfilvcf op^ii'&'ting 1ft inlnfi 
6ftvlil"Ofttt6ftt la atroftgly Iftfluanced by fton-Gaiiaaiatt atrnioapharic 
ftolae which la Impulalve In nature, and the power line harmonica. 
The power line harmonica which are integral multlplea of power 
frequency have atrong Interfering effect upto several kilohertz. 
In general, many approachea can be taken In evaluating receiver 
performance for a non-Gauaalan nolae, baaed on an appropriate 
nolae model [10]. Two approachea are uaually taken to develop a 
nolae model at ELF/VLF, depending on the bandwidth of the 
receiver. At VLF the ratio of bandwidth to centre frequency ia 
uaually amall and therefore the received nolae may be conaidered 
aa a narrow band proceaa. At ELF thla ratio ia not amall, and the 
nolae ahould be treated aa a wide band proceaa. 

Omura and Shaft [11] have conaidered VLF atmoapherie nolae aa 
a narrow band proceaa with a log-normal envelope of the form 

act) « A aln (w^t + e) 

Where n(t) la a zero-mean real atatlonary Gauaalan proceaa. The 
exponential function emphaalzea large amplitudea which have 
greater influence on the performance of a communication aystem. To 
evaluate receiver performance with thla noiae model, they 
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separately consider long and short integration times of the 
receiver. In case of long integration time the noise bursts get 
averaged out so that the central limit theorem can be applied. 
This allows one to consider the noise to be a Gaussian process. 
For short integration time intuitive as well as rigorous methods 
may be used. Though, the intuitive methods, in general, do not 
take into account the actual receiver structure, they can yield 
numerical results which are in good agreement with experimental 
values. However, It can not be Justified on theoretical grounds. 
The rigorous methods are difficult to carry out by exact 
mathematical procedures. For example, in case of a PSK receiver 
with input signal •/? alnei^t, the output is given by 

T 

V (T) » 2PT + 21 A co&(u t+d) VSP sin w t dt. 

o J o o 

0 

The evaluation of this Integral la very difficult due to the fact 
that the variables n(t) and 0 may change over the integration 
period tO,T]. To evaluate this Integral they assume T to be 
sufficiently small such that n(t) and & are essentially constant 
over this time Interval. This assumption is usually not true. 
However, they assume this in order to complete an otherwise 
extremely difficult mathematical problem. The Justification for 
this assumption la given from the point of view of end results 
only, which are reasonably in good agreement with experimental 
values. 
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Field and Lewinsteln ri?i j ^ , 

nave modelled atmospheric noise as 

a two component random process. For ELF a wide-band approach and 

£qr VLF a narrow-band approach have been taken. The basis for the 

two component representation lies on the fact that wideband 

observation of ELF and VLF atmospheric noise reveals presence of 

intermittent, non-overlapping, large amplitude spike like pulses 

superimposed on a homogeneous background. These pulses are caused 

by local thunderstorm activity while the background is caused by 

large number of relatively weak, unresolved pulses from distant 

lightning flashes. They consider the background noise to be a 

zero-mean Gaussian random process, and assume a power-Rayleigh 

distribution for the Impulsive component. The tails of the 

power-Rayleigh distribution falls off much more slowly as compared 

to a Gaussian density function. The model thus allows the 

relatively frequent occurence of large amplitude pulses. The two 

components are convolved to obtain the probability distribution 

function for the composite noise at the receiver Input. The 

numerical results for noise amplitudes are In good agreement with 

experimental values. The amplitude density function can be 

directly used for receivers which employ envelope detection. 

However, for a non-envelope type receiver e.g., a coherent PSK 

receiver the probability density function at the receiver output 

should be known for computation of bit error rate. The methodology 

to accomplish this from the Input density function has not been 

discussed' by the authors . Ih this chapter we present a method to 

compute bit error rate for PSK receiver in atmospheric-type noise. 

We also consider the effectiveness of ^a* direct .Sequence spread 
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,pe=tru» fecelvpp In preaence of atnoaphorio nolaa. alone vlth an 
Intecferlne toti®. 

B. a A NOISaE MODEL FOR ATMOSPHERIC NOISE 

Ue attempt to develop a model for the atmospheric noise which 
can reasonably describe the actual noise process, and at the same 
time allows to evaluate the receiver performance in a rigorous 
manner. The model la based on the wide band approach and It 
considers the atmospheric noise as a two component random process. 
One component la the background noise which is described by a 
zero-mean Gaussian process. The other component represents the 
non-over lapping pulses and it la considered as a discrete random 
process. These pulses although very spiky, have a finite width, 
but for convenience, these pulses are assumed to be pure impulses. 
The strength of an impulse la defined as the area under the finite 
width pulse which it replaces. The arrival of the pulses la 
assumed to be a stationary Poisson process. The pulses can be 
either positive or negative but they are assumed to be of same 
strength. This can be viewed as if the spikes are being clipped at 
a level much higher as compared to the Gaussian component. The two 
components are specified as follows. 

Background nolso 

The background noise is assumed to be additive Uhlte Gaussian 
nola* (AUGN) n^, vith » two sl<l«d apactral density of 11^/2. Thia 
InpHaa that tha Gauaaian nolaa haa a apactraa that ia auoh wldat 
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than the bandwidth of the receiver and its spectral density is 
essentially constant over the receiver bandwidth. 

Ijnptdsive component 


The impulsive component i(t) is represented by a sequence of 
impulses arriving, at random points in time and can be written as 

00 

i(t) - J act (5.1) 

n**! 

where ^SCt) is the Dirac delta function, A is the amplitude of the 

n 

impulse at time A given impulse can be positive with 

probability P and negative with probability P , i.e.. 


Prob (A^ “ 

Prob (A cn -A ) = P“ = Cl-P'^) 
n g 


(5.2) 


where A is the magnitude of the impulse. The Poisson arrival 

a 

process for the impulses can be written as: 

p(N) a exp(-\T) for N = 0, 1 , 2 , . . . . (5.3) 


where P(N) is the probability of arrival of exactly N impulses in 
Unse duration T and \ is the average arrival rate. 
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This noise model is used to evaluate 

1,0 evaluate the performance of a 

coher^ei^t PSK recelvet' In atmoapheric noise. A conventional PSK 
t-eceivei^ can not perform satisfactorily in the presence of a 
strong continuous wave interference, for example, a power line 
harmonic tone. In such situations spread spectrum or other 
communication techniques have to he employed. A direct sequence 
spread spectrum receiver with PSK signal format can be treated as 
a conventional PSK receiver to evaluate the receiver performance 
in Gaussian as well as in atmospheric noise. A direct sequence 
spread spectrum transmitter with binary PSK signal format needs 
less bandwidth than for frequency hopped scheme. Since, power 
saving la a critical requirement for a trapped miner, the former 
has a decisive advantage over the latter from this pont of view. 
The following cases are considered. 

(1) A coherent PSK receiver in absence of Interfering tone, 

(il) Direct sequence spread spectrum receiver in presence of an 
interfering tone at the carrier frequency. 

S.3 PSK RECEIVER 

Figures 5. la and 5.1b show a base-band data bit stream and 
the corresponding PSK modulated carrier wave. The base-band signal 
d(t) is polar and it takes the values of either 1 or -1 
corresponding to the binary symbols 1 and 0 respectively. The 
symbols 1 and 0 are transmitted as A^cosco^t and -AgCosto^t 
respectively. Fig. 5.1c describes a coherent PSK receiver. The 
receiver multiplies the received waveform by the synchronised 
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(Q) 





(c) 


Fig. 5.1 (a) Base band data bit stream . 
( b) PSK modulated carrier wave. 
( c) Coherent PSK receiver . 



replica of the unmodulated carrier. Ue assume a perfect 
aynchronizat ion for the carrier. 

Let, “1 be the transmitted data symbol. The received signal 
can be written as : 

Sj^Ct) = -A^ cos w^t + N^(t) (5.4) 

vhere N^(t) la a random variable which represents the atmospheric 
noise. The receiver output la given by 

T 

V^^(T) « -2A^ I cos^6i^t dt + N^(T) 

•^0 

T 

=> -A j (cos 2<o t + 1) dt + N.(T) 

S I 0 A 

^0 

= -A^T + N^(T) C5.5) 

where the double frequency component 2<.>^ has been ignored. The 
random variable Nj^(T)la due to the atmospheric noise. Probability 
density function (PDF) of N^(T) is to be known to compute the 
probability of bit error. Assuming a data symbol to be equally 
likely to be ± 1, the error probability p^^ can be written as : 

p.i “ ll'V * 

Ue now con«i4er the evaluation of the PDF of N^(T). 


> 0 


(5.6) 
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S.3.1 PDF for N^CT3 

Let Pj^C^c) and p2(y) be the PDFs at the integrator output 
due to the Gaussian and the impulsive component of atmospheric 
noise. Convolution of p^(x) with P2(y) gives the PDF for the 
random variable N^(T). One can write the resulting distribution 
function P3(z) as : 

00 

pgCz) = I P^Cz-y) p2Cy) dy (5.7) 

-00 

where, z represent the random variable Nj^(T). The PDF Pj^(x) is a 

2 

zero-mean Gaussian density function with variance & - N^T and can 

be written as: 

p (X) = -? . (5.8) 

Integrator output due to a single impulse 

Let an impulse arrives at time t^ during the integration 
period [ 0 ,T]. The integrator output due to this impulse is 
given by 

T 

V a r ± A <5(t-t^)2cos6i^t dt 

P J ® 

0 

» ± 2 A^ cosei^tj^ . 


(5.9a) 
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«e consider the worst case situation, and replace cosoi^t by a 

rectangular wave form C^(t) as shown in Fig. 5.2. With this 

replacement the output due to the impulse can be either +2A or 

g 

■■2A I depending on the polarity of the impulse and the polarity of 
the half cycle (pulse) of C^(t) during which the impulse occurs. 
Defining, “ ^^g' output due to a single impulse can be 

written as 


''p ‘ 


(5.9b3 


The following possibilities can happen. 


(a) The impulse la positive and occurs during a positive pulse: 

V • 

(b) The impulse is negative and occurs during a negative pulse: 

V « A« . 

P G 

(c) The impulse is positive and occurs during a negative pulse: 

V =» -A . 

p 


(d) The impulse la negative and occurs during a positive pulse: 


« -A^ 
P G 


Assuming there are equal number of positive and negative half 

cycles in the wave form C (t) In the interval [0,T], it la equally 

s 

likely that the Impulse will occur during a positive or negative 
keif cycle of C^(t). Given that an impulse has arrived, one can 

deduce the following 
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Fig. 5.2 (a) Local sinusoidal carrier waveform 2cos^ot 

( b) Rectangular waveform Cs(t) used in place of the 
local carrier sinusoid for computational purposes. 



Prob[v^ = A,] = I P* » 1 p- 

1 

"■ I ' (5.10) 

Prob[Vp = -aJ = ‘ p- * 1 p* 

1 

“ 1 • 

From (5.10), It la aeen that the output for a given Inipulae la 
equally likely to be ±Aq, and it la independent of any probability 
that may be attached to the impulae being poaitive or negative. 
Thua it allova ua to aaaume that all the iapulsea are of positive 
polarity without any loss of generality. Uhen there are more than 
one impulae (during the time interval [0,T]), this argument atlll 
holds . 

Integrator output due to any number of impulses 

Suppose two impulaea have occured during the interval [0,T]; 
the output can be 2 Aq, 0 and -2Ag, which can be Illustrated by the 
following poaslble combinations : 
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The flfst term Inside the brackets is the output due to the first 
luipuls* and the second one is due to the second impulse. Each of 
the above combinations takes place with probability -j. In a 
similar way output for the arrival of three impulses can be 
obtained. When there are exactly one, two or three impulses, the 
conditional probabilities can be written as follows*. 

Cal Number of Impulses ■ 1 


Output 


0 

-^G 

Prob 

PCD 

2 

0 

PCD 

2 


Cbl Number of Impulses ■ & 


Output V 

P 

2^0 

A G 

0 

-^G 


Prob 

PCD 

i 

0 

PC2) 

2 

0 

PC2) 

4 


Ccl Number of impulses ■ 3 


Output 

3Aq 


^G 

0 

-^6 



Prob 

PC3? 

8 

0 

3PC3) 

8 

0 

3PC3) 

8 

0 

PC3) 

8 
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where P(1), P(H) and P(3) are probabilities 
one, two and three Impulaea during the 
followine observatlona can be made from thes 


of arrival 
interval 
e examples 


of eicactly 
[0,T]. The 


(») Prob (Vp = KAg) = P„b (Vp = -KAg), vhe<-* K = 0,1, 2, 3,... 

(b) Pop' N nurobec o£ Impulaea the fliaxiniuin value of K is equal to N 
or K ^ N. 


(c) 


When N ia even 
> 0, occur for 


(odd), non-zero probabilities 
even (odd) values of K. 


l.e.ProbCV =KA_) 
p G"^ 


For inteerator output to be equal to KA^, when N impulses 
have arrived, the number of Impulses each contributing A^ is equal 
to where impulses contribute -Ag (by each of them), 

aubject to the even-even and odd-odd relationship between K and N. 
For example to get output equal to Ag (K=l) when there are 3 
impulses, one can have the following combinations, where in each 
of the triplet there are two Ag and one -Ag. 

|(Ag,Ag,-Ag) AgJ, {(AQ,-Ag,Ag) Ag}, |(-Ag,Ag,Ag) + Ag} . 

Ih general one can write the conditional probability for = KAg 
when there are N Impulaes, as follows. 

lotting u « 
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ProbCKAgIN) 


Prob {" => C'-K) => -v4 

^ 

Prob [N arrivals] 


where u => +V6 meana u number of impulses occur during the period 
of positive half cycles of C^(t). Similar is the meaning for (N-u) 

a 

s> -ve. Hence, one can write 


Prob (KAqIN) 


6Xp(“AT/2)(XT/2)‘* 6xp(-XT/2)CXT/2)“ “ 

u i (H-u) i 

exp(-XT) (XT)“/N ! 


SB 


S3 



(5.11) 


Hence, the unconditional probability for = KA^ (with K ^ N) is 
given by 


oe 

Prob (KAq) = ^ P(N) 

N=K+2i 


N 

N+K 

2 


(i)“ 


(5.12) 


0 , 1 , 2 , 3 . - . 


This la a discrete probability distribution function with the 
symmetry Prob(KAQ) » Prob(-KAg)- The PDF defined by (5.12) was 
earlier denoted by section (5.3). 
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Comparing (5.7) with (5.12). the distribution function P 3 (z) 
for the random variable N^(T) at the integrator output due to 
atmoapheric nolae can be written aa : 


p3(z) 


OO <0 




Pj(z-6(y-KAQ)) P(N) 


N 

K+N 

2 J 


(j) dy 


(5.13) 


i - 0,1,3 


where Pj^(3C) and P(N) are given by (5,8) ai^d (5.3) respectively. 
The integral (5.13) for a given value of K is equivalent to 
shifting the Gaussian PDF given by Pj(x) to the position KAg and 
multiply by the factor fj^ which can be seen from (5.13) to be 



The PDF P^la) la an infinite summation of such weighted and 
shifted Gaussian shaped PDFs. The bit error rate (BER) for a PSK 
receiver can be computed by using (5.6). Before computing the 
numerical results for a PSK receiver, let us formulate for the 
error probability in a direct sequence spread spectrum receiver. 

S.4, DIRECT SEflUENCE SPREAD SPECTRUM COMMUNICATION SYSTEM 

Fl*. B.3 ohowa a direct aequence (OS) apread apactrum 
coaauplcatlon ayataa, with binary PSK aa the nodulatlon fornat. In 
the trananltter the carrier containing data atrean la nultlpUed 
hy a pa.udo random or p.eado nol.e <FN) aequence. The PN a.quence 
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changes ita atate many tlmea during a data pulse o£ duration T. 
This leada to piiaae shifting o£ the carrier at a snuch higher rate 
than that which would have been caused by a data stream alone. The 
multiplication process by the PN sequence can be Identified as to 
spread the spectrum of the signal. The receiver multiplies the 
received signal with a synchronized replica of the transmitter's 
PN sequence. This remultiplication at the receiver despreads the 
spectrum of the signal which now appears in the original bandwidth 
of the data stream. The subsequent recovery of data Is identical 
to that of a conventional PSK receiver. 

Figures 5.4a and 5.4b describe the relationship between a 
data bit and the PN sequence which is being generated from a four 
stage maximal length shift register shown in Fig. 5. 4c. If the 
register has N stages, one can write the following: 


(1) 

Length 

of the sequence = 

2“ - 1. 


(11) 

No. 

of 

ones in the sequence = 

2<N-1) . 

(5.15) 

(ill) 

No. 

of 

zeros in the sequence = 

1 

1 

s; 



E»ch titat* o( th« PN la called a chip and the chip 

deration i. danot.d by I„. H a data bit la I a.cond long, than 

T 

T a — seconds. 
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Flg.5.3 Adirect scq.uence spread spectrum 
communication system- 








Ill 




If at t*0;thc register content is 0001 (left to right)the sequence is 
10001111010110O 

t t 


Fig.5.4 (a) Data bit. (b)PN sequence. (c)PN sequence 
generator r 
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S.4. 1 Error probability : 

Let, ~1 be the transmitted data symbol. The received alfinal 
Is given by 

S2Ct) = -A^ PN(t) coa<o^t + N^(t) + ICt) (5.16) 

where la the amplitude of the signal and PN(t) la the pseudo 
noise aetjuence. The variables N^(T) and I(t) represent atmospheric 
noise and the power harmonic Interfering tone respectively. The 
interfering signal with amplitude a. and frequency can be 

written as 

ICt) = Oj coa(o3j^t + Sj) (5.17) 

where dj Is the phase angle of the Interfering tone. The phase 

angle <9j Is a random variable which la assumed to be uniformly 

2 2 

distributed over the range of [0, 2n]. The ratio "j/A^ indicates 
the strength of the Interfering power In comparison to the signal 
power. In decibels (d(3) this ratio which la usually denoted by 
J/S, is given by 

s - 

s 

Asauming a perfect synchronization for the carrier and the PN 
sequence, the receiver output written as 

VajCT) - -A.T 
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where -A^T^Nj^CT) and fi|(T) are the outputs due to the data symbol, 
the atmospheric noise and the Interfering tone, respectively. The 
PDF for N^(T) Is still given by (5.13), which has been obtained 
for a PSK receiver. In the present case, the multiplication of the 
received signal by the PN sequence results In phase shift keying 
during a data symbol , but the basic requirement that there must 
be equal nuaiber of positive and negative half cycles of the local 
carrier during the Interval [0,Tj remains the same as that of a 
PSK receiver. For a given phase angle d, the variable g.(T) is 
given by 

T 

gjj^(T) “ J PN(t) co8|(<o^-<o^)t + dt 

» B slnd + « B cosfif (5.20a) 

X 0 1C 


where 


B^ - J PN(t) sin(<aj“ai^)t 


(5.20b) 


4. 

« r PN(t) cos(6)^-<o^)t . 
■^0 


(5.20c) 


The variable Sj^CT) depends on 


the PN code used and the position of 


th. int.rf.rln* ton. with r..p.ct to th. c.rrl.r, U. con.ider th. 


case of « m and hence one can write 
1 o 


B 


« 0; « J BN(t) dt 

a 0 


(5.21) 
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From (5.15)y (S.SOa) - (5.S1), we get 

=> otj cose . (5.22) 

The receiver output for ej. = e la eiveti by 

^a2^^^ = -AsT+ N-CT) + a^cosei (5.23) 

A 1 C 

e 

can be 


(5.24a) 


(5.24b) 


8.B NUMERICAL RESULTS 

For PSK receiver the bit error probability p^j^ for a given 

signal to noise ratio (SNR = '^b'^^o’ ^b ~ computed 

in terras of the pararaeters \T and a new parameter g^ which is 
defined by the relation AJ la the receiver output due 

to a data bit while is the output due to a single impulse. Thus 
g^ is a measure of the impulse strength in terms of the output due 
to a data bit. The parameter XT can be Identified as the average 
number of impulses arriving during the Interval [0,T]. 


Hence, the conditional probability of error for d , 
written as 


Probj^p^2j^j “ Pit'ob - '“AgT+ N^(T) + oijCos^^ T^ ■ > 0 


The unconditional probability of error ia given by 


Prob 


[^62 



Figures 5.5 to 5.7 give p 
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for different values of XT and 
The curve with = 0.0 la £or the case when there la no 

impulsive component. For b^= 0.1 and XT 5 0.5, the additional error 
rate due to the Impulsive component is quite small. The 

performance degradation is within 0.5 dB of SNR. Uith this value 
of if XT=1.0, the performance degradiation is within 1.5 dB. 
The amplitude g^ has much greater effect on the receiver 

performance than the parameter XT. This is because, if the number 
of pulses arriving at the receiver input increases, the output due 
to the pulses does not increase proportionally; there are 
cancellation effect. On the other hand an increase in g (XT 

X 

fixed) does not lead to any further Increase in the existing 

cancellation level. This can be seen by comparing different curves 
with same gj|.^T product. This fact demonstrates that for a given 
impulse strength and arrival rate, the effect of the impulsive 
component on the receiver performance, can be reduced by 
increasing the integration time. 

Fig. 5.8 describes the error probability for a direct 

sequence spread spectrum receiver in presence of an interfering 
tone at the carrier frequency. The PN code is derived from a 
maximal length shift register of stages, resulting in a sequence 
length of 63. The receiver performance neither enhances nor 
deteriorates in presence of the atmospheric noise when there is no 
interfering signal; its performance la same as that of a PSK 
receiver. However, it provides strong immunity to the interfering 

aignal. For example for J/S upto 10 dB, the peformanee of the 
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receiver due to the intefering tone is marginally affected; the 

error degradation is within 0.2 dB of SNR. For J/S = 20 dB and p 
5 62 
10 , this figure is within 1.5 dB. 

S. 6 A DISCUSSION ON THE ASSUMPTIONS MADE REGARDING NOISE MODEL 

Two appr oximat i ons have been made in this chapter for 
deriving the PDF for the Impulsive component. These approximations 
are pessimist Ic in nature for the following reasons. 

(i) The finite width pulses have been replaced by pure impulses. 
In the case of finite width pulses there will be cancellation 
effect by the oppositely polarized neighbouring half cycles 
of the local carrier. Because some part of the pulse may be 
on a positive half cycle and the remaining part on a negative 
half cycle, this will result in a reduced integrator output 
for the impulsive component. 

(li) The sinusoidal local carrier has been replaced by a 
rectangular wave form C^<t) to obtain the receiver output for 
the impulsive component in Eq. (5.9a). Since, the amplitude 
and frequency of C^C^:) are equal to that of the original 
sinusoid, the receiver output values (in magnitude) are 
always greater or equal to those values which would have 
been obtained by the use of the pure sinusoid. 



Probability of error 
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Fig.5.5 Performance of the PSK receiver in atmospheric 
noise 



Probabil'ity of error 
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Fig ,5.6 Performance of the PSK receiver in atmospheric 
noise. 




Probability 



c‘ c »7 o ^ the PSK receiver In atmospheric 

Fig. 5.7 Performance ot r^rv 

noise. 




Probability of error 


120 



Fig. 5.8 Performance of the spread spectrum receiver 
in atmospheric noise and tone interference 



CHAPTER 6 


CONCLUSIONS 

A loop antenna located on the earth aurface and operating at 
VLF/ELF banda , uaed for underground mine communication , haa been 
analyaed by the method of momenta. The analyaia doea not require 
any numerical Integration aa all Integratlona aaaoclated with the 
solution procedure have been expreaaed in cloaed forma. It has 
been found that aa the loop diameter la increased (keeping other 
parameters fixed) behaviour of the loop l.e., variation of input 
admittance and current distribution change very rapidly when the 
magnitude of the input auaceptance is greater than about seven 
times the value of the input conductance. This means, the 
transition from a small loop behaviour to that of a large loop la 
very rapid. Uith the change in current diatribut Iona , the field 
inside the earth also changes accordingly. For a small loop the 
vertical magnetic field component is almost independent of the 
relative angular position of the observation point with respect to 
the feed point of the antenna. For a large loop the magnitude of 
the fields depend both on radial distance from the loop axis and 
the angular position of the observation point. Uhen there is a 
large variation in current distribution the radial component along 
the line Joining the feed point and the loop axis is very large 
compared to radial or a vertical components of the magnetic field 
at any other angular positions. It is useful to predict how large 
a loop can be treated to be a small loop without going through the 
method of momenta. A simple 


expression for this 


has been 
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developed. It is to be noted that it is possible to compute the 
input resistance of a small loop in free space by assuming that 
the loop carries an uniform current. However, it is not true for a 
loop antenna on the earth surface or one embedded in a dissipative 
medium. The reactive part of the input impedance for a small loop 
can be assumed to be same whether the loop la in free space or 
located on the earth surface. 

A model for atmospheric noise which contains an impulsive 
component; superimposed on a Gaussian background has been 
presented. Performance evaluation for the conventional as well' as 
spread spectrum receiver with PSK signal format sub Jected to 
this noise has been dlseuased. It has been shown that same 
procedure is applicable for both conventional as well as spread 
spectrum receiver. It has been found that when the area under a 
pulse (impulse strength) associated with the impulsive component 
is less than one tenth of the area under a data bit, the 
performance degradation is within 0.5 dB of signal to noise 
ratio. This figure la baaed on the assumption that average arrival 
rate of impulses is ^ 0.5 during the period of a data bit. 


In the present Investigation we have considered the earth as 
a homogeneous lossy half apace. A better model will be to treat it 
as a stratified media. For a stratified medium it may not be 
feasible to express the potentials in a closed form. However, the 
potentials can be expanded in a series form by taking some 
approximations. There after the analysis of the antenna problem 
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can be carried out on similar lines as described in this investi- 
gation. Similarly, the self Impedance for a amall loop antenna 
over a atratifled earth can be obtained. 

In mine noiae environment due to preaence of strong 
interfering signals spread spectrum communication technique will 
be useful. Spread spectrum communication system can also be used 
for multiaccess capability using code division multiplexing. The 
performance analysis of the spread spectrum communication system 
in VLF/ELF bands will be along the usual lines as for high 
frequency environment, with the noise model developed in the 
present Investigation being used in place of Gaussian noise. 
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